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Abstract 

We propose a dual formulation for the S Matrix of M = 4 SYM. The dual provides 
a basis for the "leading singularities" of scattering amplitudes to all orders in pertur- 
bation theory, which are sharply defined, IR safe data that uniquely determine the full 
amplitudes at tree level and 1-loop, and are conjectured to do so at all loop orders. 
The scattering amplitude for n particles in the sector with k negative helicity gluons is 
associated with a simple integral over the space of k planes in n dimensions, with the 
action of parity and cyclic symmetries manifest. The residues of the integrand compute 
a basis for the leading singularities. A given leading singularity is associated with a 
particular choice of integration contour, which we explicitly identify at tree level and 
1-loop for all NMHV amplitudes as well as the 8 particle N 2 MHV amplitude. We also 
identify a number of 2-loop leading singularities for up to 8 particles. There are a large 
number of relations among residues which follow from the multi-variable generalization 
of Cauchy's theorem known as the "global residue theorem". These relations imply 
highly non-trivial identities guaranteeing the equivalence of many different representa- 
tions of the same amplitude. They also enforce the cancellation of non-local poles as 
well as consistent infrared structure at loop level. Our conjecture connects the physics of 
scattering amplitudes to a particular subvariety in a Grassmannian; space-time locality 
is reflected in the topological properties of this space. 



1 A Dual Theory for the S Matrix 



Scattering amplitudes in gauge theories and gravity possess extraordinary structures that 
are completely invisible in the textbook formulation of quantum field theory. Much of the 
recent progress in understanding scattering amplitudes was triggered by Witten's twistor 
string theory for Af = 4 SYM amplitudes [TJ. The stimulus provided by this proposal 
led to a number of powerful new methods for explicitly computing amplitudes, many of 
which arose from a fresh examination of older ideas in field theory and S Matrix theory. 
These include the BCFW recursion relations for tree amplitudes [21 [3j HI [5], El El [8] , along 
with their supersymmetric extension [9], [10] . Amongst other things, the recursion relations 
and their explicit solution |1TJ have allowed for a direct verification of the remarkable dual 
conformal symmetry [T2l H~3] of SYM amplitudes at tree level [14.]. Furthermore, as was 
realized quite early on by Hodges [T51 [TBI HT] and rigorously understood more recently in (2,2) 
signature [HI [19], the BCFW recursion relations are most naturally formulated in twistor 
space [20, EH E2J E3], where a beautiful "Hodges diagram" formalism exposes structures 
in the amplitudes that are obscured in momentum space. All of these amazing properties 
strongly suggest the existence of a "weak-weak" dual theory which directly computes on-shell 
scattering amplitudes without resorting to the conventional evolution through space-time. In 
this paper we propose a dual formulation for the S Matrix of = 4 SYM along these lines. 

1.1 What Should the Dual Compute? 

As a preface to our proposal, it is worth discussing what properties should be expected of a 
dual theory for the S Matrix on general grounds. The most basic question of all is simple: 
what should such a theory compute? At tree level, scattering amplitudes are well-defined 
rational functions of the kinematical invariants, but loop level amplitudes suffer from infrared 
divergences. Is there any sharply defined data about loop amplitudes that might serve as the 
output of a putative dual theory? 

A look at the general structure of scattering amplitudes at 1-loop offers a very natural 
answer to this question. As is well-known, any 1-loop amplitude in four dimensions can be 
written as a linear combination of a basis of scalar integrals with rational coefficients [21]. A 
generic theory will have scalar box, triangle and bubble integrals, as well as purely rational 
terms. All of the infrared divergences are isolated in the basis of scalar integrals; the rational 
coefficients then carry all the information about scattering amplitudes at 1-loop, and can 
therefore serve as sharp "data" to be computed by a well-defined theory. 

As reviewed and discussed at length in [9], these rational functions have a beautiful 
and very physical interpretation. Loop level amplitudes have branch cuts as a function of 
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the kinematical invariants, and so it is natural to examine the discontinuities across these 
branch cuts. The discontinuities themselves can have branch cuts as a function of other 
kinematical invariants and the procedure can be continued, finally arriving at an object 
which is aptly named the discontinuity across a leading singularity, which is the highest co- 
dimension singularity possible. In what follows we will often abuse nomenclature and use 
the phrase "leading singularity" to mean "discontinuity across the leading singularity" . Just 
as the textbook unitarity cut is computed by putting two propagators on-shell, the leading 
singularity at 1-loop in four dimensions is computed by putting four propagators on-shell 
[25]. 

A related point of view is to consider the sum of Feynman diagrams in a 1-loop amplitude 
which share four given propagators. For an n-point amplitude, the four propagators in the 
loop separate the n-particles into four sets. Let K~x, K 2 , K 3 and K± be the sum of the momenta 
in each set. It is possible to put all these Feynman diagrams together as a single integral of 
the form 



where we have exhibited the four special propagators and where R(£) is a rational function 
of £ and the external momenta which contains the rest of propagators and tensor structure 
in the Feynman diagrams. The integral ([1]) has IR divergences and in a general theory 
it could have UV divergences as well. All divergences come from choosing the region of 
integration to be M 4 . If, instead, we choose to interpret the integral as a contour integral 
in C 4 with contour T 4 - (S 1 ) 4 given by \£ 2 \ = e, \{£ - i^) 2 | = e, \(£ - K x - K 2 f \ = e and 
\(£ + K 4 ) 2 \ = e with e a small positive real number, then the integral becomes finite. In 
mathematical terms, the new contour integral is computing a residue! In order to compute 
the residue we have to evaluate R(£) on the points, £*, in C 4 which satisfy the equations with 
e = 0. There are two solutions and therefore two different T 4 contours. The two contours are 
related by parity. Physically, R(£*) factors as the product of four tree amplitudes. Explicitly, 
R(£*) — > M^] e e M^ ee M^ ree M^ ee , where the subindices mean u =upper, I =left, r=right, 
and b =bottom, to refer to the location in the figure below. The sum is over all physical 
states in the theory which propagate in the internal legs [25J. 




(1) 
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R(£*) is a leading singularity of the amplitude. At 1-loop something special happens. In 
the basis of scalar integrals, the scalar box integrals are the only objects that can be cut 
in this way (with suitable normalization giving "1"), and so the rational coefficients of the 
scalar box integrals are precisely the leading singularities at 1-loop [25]. More precisely, the 
coefficient of a given box integral is the sum over the two leading singularities associated with 
the two different solutions for £* or equivalently the two T 4 contours. 

An ansatz for the amplitude which includes only scalar boxes with these coefficients is 
an excellent initial guess for the whole amplitude, since it is guaranteed to reproduce the 
leading singularity. This ansatz also yields contributions to subleading singularities, which 
involve cutting fewer lines. In a generic theory, the set of subleading singularities produced 
in this way is incomplete, and therefore the lower-point scalar integrals and rational pieces 
must be present to compensate for the difference. The remarkable feature of maximally 
supersymmetric theories at 1-loop is that fixing the correct leading singularity suffices to fix 
all the subleading singularities as well, so that only the scalar box integrals appear [26J: 






The procedure for moving on to higher loops is now clear [271 |28| [29] . At /-loop order one 
can choose Feynman diagrams with a given set of 4/-propagators, with each loop variable 
appearing in at least four of the propagators of the set. The loop integral 

d% . . . d% 



Ylf =1 P^£i,...,ii,h,...k n ) 



i- 



(2) 



which generically suffers from IR divergences (UV divergences are absent in M = 4 SYM), 
becomes a completely well defined object if we take it to be a contour integral in C 4 ' with 
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a T Al contour given by |P X 2 | = ... = |P|J = e. Once again, R(£i, ...,£/) is the product of 
tree amplitudes summed over all possible physical on-shell states in the theory propagating 
in the internal lines. As a simple illustration consider the following 2-loop 5-point leading 
singularity 




where the open circles denote tree amplitudes and all internal lines are on-shell. 

A comment is in order here. For a given number of particles, n, it is impossible to find 
diagrams with 41 propagators if I > n — 3. This problem is resolved by defining "composite" 
leading singularities. The theory of composite leading singularities was developed in [27J ES 
l30l 151] . At the end of section 3.2 we illustrate it with a four-particle amplitude at 2-loop 
order. 

Leading singularities, including composite ones, constitute perfectly well-defined data 
associated with scattering amplitudes at all orders in perturbation theory. It is natural to 
conjecture that, for maximally supersymmetric theories, the leading singularities continue 
to fully determine the amplitude at any loop order as well; for M = 4 SYM, some strong 
evidence for this proposal has been provided at 2- and 3- loop order [27J, [30l EU [32] . We will 
optimistically assume that this conjecture is true, and that the leading singularities not only 
constitute well-defined data about the S Matrix at all orders in perturbation theory, but also 
suffice to completely determine it. 

It is important to mention that the direct connection between leading singularities and 
scalar integral coefficients at one-loop does not generalize to higher loops. At higher loops, 
there is not a unique basis of integrals. In the literature of higher loop M = 4 amplitudes 
the main emphasis has been on expressing loop amplitudes as a sum over scalar integrals 
and generalizations, which include loop momentum dependent numerators, with coefficients 
which are rational functions of the kinematical invariants (see e.g. [341 [29]). In [30], it was 
proposed that the coefficients appearing in the linear combination of scalar integrals can be 
completely determined by correctly matching all the leading singularities. This has been 
tested for 5 and 6-point amplitudes at 2-loop and 5-points at 3-loops [30j EU I52"J*1 . It is not 
the purpose of this paper to determine amplitudes in terms of a given basis of scalar integrals; 

lr The five-particle and the parity even part of the six particle 2-loop amplitudes where originally computed 
using the unitarity-based method in [33] . 
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rather we aim to determine the leading singularities directly, which are physical and basis 
independent. Only at one-loop will we interchangeably talk about leading singularities and 
scalar boxes, with the hope that this will not generate confusions when discussing two-loop 
leading singularities. 

Summarizing, we have arrived at a natural picture for what a dual theory for the S Matrix 
should compute: it should produce the leading singularities of scattering amplitudes at any 
loop level. 

1.2 What Symmetries Should It Manifest? 



We will use spinor-helicity variables A a , A a associated with the momentum p a as p" a = A" A" 
where a — 1, • • • ,n indexes the n external particles [35]. Since all helicity states are in the 
same supermultiplet, we can also label the external particles by Grassmann variables fj a so 
that the amplitude is a smooth function of (A a , A a ,f/ a ). Now, the usual local formulation of 
M = 4 SYM makes two symmetries manifest. The Bose symmetry of the non-supersymmetric 
theory is enhanced to a cyclic symmetry of the color-stripped superamplitude [36], which 
simply acts as g : a — > a + 1 (mod n). The other symmetry is parity. 

Nonetheless, even at tree-level both of these symmetries are obscured by the BCFW 
formalism. Manifest cyclicity and parity are broken due to the choice of two special external 
lines and the deformation of A^ and Xj t while leaving untouched. Since cyclic 

symmetry and parity are not manifest, they appear as highly non-trivial identities relating 
different BCFW representations of the same amplitude. 

The simplest non- vanishing amplitudes in SYM theory are in the sector with two negative 
or positive helicity gluons, the so-called MHV and MHV amplitudes. In this case the super- 
amplitude is determined by the famous Parke- Taylor amplitude [2Z] , and is manifestly cyclic 
and parity symmetric. The first example illustrating the tension with manifest cyclic and 
parity symmetries occurs in the sector with 3 gluons of a given helicity (the NMHV ampli- 
tudes) and the simplest case is the 6 particle amplitude. Consider for instance the alternating 
helicity configuration. Applying the usual BCFW deformation on Ai, A2 yields a three-term 
representation for this amplitude, which manifests the (obvious) symmetry g 2 : a — > a + 2 
shifting the particle index by two units: 



^ CF + w + - = (1 + 9 2 + 9*" 



(46) 4 [13] 4 



L[12][23](45)(56)(p4 + Ps + Pe) 2 (6|5 + 4|3](4|5 + 6|1]J 



(3) 



Now with the supersymmetric form of the BCFW recursion relations, we can also consider 
the parity conjugate of this deformation, where A2, Ai are deformed. This yields a parity- 
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conjugated form of the same amplitude given by 



^p + ( b£fw) = (1 + 9 2 + g 4 



[3|(2 + 4)|6> 4 



(4) 



[23] [34] (56) (61) (p 5 + p 6 + Pi) 2 (1 16 + 5|4] (5|6 + 1 12] J 

This is not manifestly the same expression as Mbcfw- The statement of parity invariance is 
instead a highly non-trivial 6-term identity 

A^BCFW — ^P(BCFW) (5) 

It is straightforward to see that the same identity is enforced by demanding cyclic symmetry of 
the superamplitude-indeed it is one component of the supersymmetric expression of cyclicity. 
At higher points, increasingly complicated identities must be satisfied: a 12 term identity for 
the 7 point NMHV amplitude, a 40 term identity for the 8-pt iV 2 MHV amplitudes and so 
on. 

The simplicity of the BCFW construction and its clear connection to deeper structures in 
twistor space suggest that this representation of the amplitudes is secretly being computed by 
a dual theory. If this is the case, it makes no sense for any particular BCFW representation 
to be privileged; we should expect that parity and cyclic symmetries are manifest in the dual 
theory. However physical amplitudes should not be output in a unique form, but should 
instead somehow be associated with equivalence classes, any given representative of which is 
not manifestly symmetric. Some powerful mathematical structure must therefore be at work 
guaranteeing the equivalence of these different representations and producing, for instance, 
the remarkable 6,12,40 and higher term identities we mentioned above. 



1.3 What Miracles Should It Perform? 

There is an even more basic reason to suspect a powerful structure at work: any dual theory 
for the S Matrix, making no reference to space-time, must provide an understanding of how 
local space-time physics emerges. The AdS/CFT correspondence [2B] has made us accus- 
tomed to the holographic generation of spatial dimensions, but the details of this mapping 
are shrouded in the mysteries of the strongly coupled CFT. On the other hand, in a "weak- 
weak" dual theory for the S Matrix, the mechanism allowing the emergence of space-time 
has nowhere to hide, and must be visible perturbatively. 

It is important to first understand how the presence of a bulk space-time with local inter- 
actions is encoded in the S Matrix. An important consequence of space-time locality is that 
the amplitude has a highly constrained singularity structure. At tree level, the amplitudes 
can only have "physical poles" , of the form 



(Pai + ■ ■ ■+Pa m ) 2 
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This is clear from Feynman diagrams. However, one can easily imagine other types of poles, 
for instance of the form 

1 (7) 
(1|2 + 3|4] y) 

These cannot possibly arise from Feynman diagrams, and do not admit any sort of local 

space-time interpretation. The most basic indication of space-time locality is therefore that 

these sorts of non-local "unphysical poles" must be absent from the amplitudes. 

Note however that the individual terms in the BCFW form of the amplitudes in equation 
([3]) do have unphysical poles which we highlighted by separating from the physical poles as 
an explicit factor. The presence of unphysical poles in each term is in perfect accord with the 
idea that these expressions are being produced by a dual theory with no regard for manifest 
space-time locality. But the particular combination of objects appearing in the full amplitude 
must ensure that the unphysical poles cancel, rendering them "spurious". Simply looking at 
equation ([3]), the cancellation of the unphysical poles looks somewhat miraculous. However, 
this cancellation follows directly from the remarkable 6-term identity M B cfw — ^p(bcfw) 
we just discussed in the previous subsection. This is because the unphysical poles appearing 
in Mbcfw are different than the ones appearing in Mp(bcfw)- 

We therefore see that the consistency of the BCFW rules with local space-time physics is 
highly non-trivial, since the individual pieces in a BCFW expression for scattering amplitudes 
can not arise from local physics. However, in this example the same powerful mathemati- 
cal structure guaranteeing the equivalence of different representations of the amplitude also 
ensures the absence of unphysical poles and makes a local spacetime interpretation possible. 

Recently, Hodges has given a remarkable interpretation of NMHV amplitudes as volumes 
of polytopes in the twistor space associated with the space on which dual conformal transfor- 
mations act as conformal transformations [39]. His representation keeps the cyclic symmetries 
manifest and beautifully explains the cancellation of spurious poles. We are seeking a sim- 
ilar understanding which does not make heavy use of the dual space, since amongst other 
things, we are ultimately interested in describing gravity, where this dual space is unlikely to 
play a central role. See also [ID] for another interesting discussion of spurious poles in tree 
amplitudes. 

Beyond tree-level, there are many other relations which must be satisfied by scattering 
amplitudes in order to be consistent with a local space-time description. At 1-loop, the 
leading double-logarithmic IR divergence of the amplitude must be proportional to the tree 
amplitude. Working in dimensional regularization, this implies |¥T1 153] 

n 

Mr ioop i IR = - 72 £(- s v+o e < ree - (8) 

i=l 
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This "IR equation" guarantees that the IR divergence of the 1-loop amplitude has the correct 
physical interpretation of accounting for the soft emission of gluons in the tree-level process. 
As we reviewed in the previous subsection, one-loop amplitudes can be written as linear 
combinations of scalar box integrals; each integral can have IR divergences that depend 
on the kinematical invariants of the particular box under consideration, and thus there are 
many IR equations associated with the double logarithms of different kinematical invariants. 
Equation (jSj) thus becomes a set of relations constraining the various scalar box coefficients. 

These equations naturally split into two types. One set constrains a combination of scalar 
boxes with double logarithms of invariants of the form (pi + Pi+i) 2 to be proportional to the 
tree amplitude, while the remaining ones constrain combinations of scalar box coefficients to 
vanish. A particular linear combination of these equations [12| [9] directly yields the BCFW 
recursion relations in their parity-symmetric form ~ [BCFW + P(BCFW)] [2] (indeed the 
BCFW recursion relations were first discovered by conjecturing that the two pieces were equal 
[2]). This has a pleasing implication. We motivated our discussion of leading singularities by 
the search for IR safe data at loop level, but of course the tree amplitudes are also perfectly 
well-defined rational functions, and it may have seemed odd to have two different kinds of 
sharp data associated with the S Matrix. Fortunately this is not the case-a subset of the 
leading singularities contain all the information in the tree amplitudes as well. 

The IR equations are a reflection of both locality and unitarity in the scattering ampli- 
tudes, and given that they directly translate into a statement about relations among leading 
singularities, a dual theory that computes leading singularities should provide a direct un- 
derstanding of them. 

1.4 Our Conjecture 

In a series of papers, we have been pursuing clues to a dual theory for the S Matrix, largely by 
studying tree- level amplitudes. In our previous paper [18] (see also [19]) we showed that the 
BCFW construction of tree amplitudes is most naturally formulated in twistor space. While 
the twistor space interpretation has been extremely useful for organizing and identifying 
hidden patterns in the BCFW representation of tree amplitudes, twistor space alone does not 
provide a formulation with manifest cyclic and parity symmetry that illuminates the crucial 
cancellation of unphysical poles. That said, in twistor space and momentum space, we saw 
that amplitudes can be written in a novel form which we called the "link representation", 
offering a hope for unifying the different BCFW terms into a single object. Pursuing this 
line of thought, we were led to consider the following object as an ansatz for the n-particle 
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tree amplitude in the sector with k negative helicity gluons: 

(12 • • • k) (23 • • • (k + 1) f • • (nl ■ ■ ■ (k - 1) ) JJ ^(^a W ( 9 ) 

where the W a are twistor variables, obtained by fourier-transforming with respect to the A a , 
W = (W\n) = (ii,\\f}),azd 

K- m fc ) = e ai - afc C aimi • • • C akmk (10) 

This expression can be thought of as a unified form of the link representation, as we discuss 
in further detail in the appendix. The first striking aspect of this formula is its simplicity, as 
the integrand is comprised of a single term. Furthermore the cyclic symmetry of the formula 
is manifest; as we will see with a tiny bit of work, parity is manifest as well. 

We will study this object at great length for the rest of this paper. Going back to 
momentum space is trivial, and the resulting expression is computed as a multidimensional 
contour integral. As usual with contour integrals, there is really no integral to be done, 
and we are instead interested in the residues of the integrand at the location of each of its 
singularities. Our tree-level motivations had led us to expect that the residues of this object 
would compute terms in the BCFW expansion of tree amplitudes, and this was correct for 
MHV amplitudes and the 6 particle NMHV amplitude. Continuing to 7 points, we identified 
all terms appearing in the tree amplitudes amongst the residues, but were puzzled to find new 
objects which were not needed for the tree amplitudes. To our great surprise, we found that 
these objects could be identified as the leading singularities of the 7 point amplitude at 1-loop! 
This also resolved a little mystery at 6 particles: the tree amplitude only involves a particular 
linear combination of residues, which left us wondering why some information was evidently 
"wasted". The 1-loop interpretation allowed us to interpret all of these residues as leading 
singularities. Moreover, while at MHV level our formula only produced the tree amplitude, 
all the 1-loop leading singularities of the 1-loop MHV amplitudes are proportional to the 
MHV tree amplitude and are thus consistent with a 1-loop interpretation. Subsequently, we 
identified residues corresponding to the 1-loop leading singularities for all NMHV amplitudes 
and the 8 particle N 2 MHV amplitude, including the coefficient of the completely IR finite 
"four-mass" boxes. 

Further investigations found the pattern repeating itself: beginning with the 8 particle 
NMHV amplitude, some of the residues did not appear in any of the 1-loop leading singu- 
larities. We were able to identify these as certain 2-loop leading singularities! Since full 2 
loop amplitudes for beyond 5 particles are not yet known, we cannot yet perform a simple 
check that all two-loop leading singularities are present, but given the striking pattern we 
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have seen it is a very plausible conjecture, which if true further suggests that the residues are 
computing the leading singularities at all loop level. One immediate multi-loop prediction 
of this conjecture is that leading singularities of MHV loop amplitudes should always be 
proportional to the tree amplitude; this is true for all cases that have been computed so far. 
Another prediction is that the leading singularities occurring in the 6 and 7 particle ampli- 
tudes at two loops and beyond should be the same as the objects occurring in the 1-loop 
leading singularities; this also appears to be the case for the 2-loop 6 particle amplitude [43] . 

It appears that we have found a rather remarkable object that knows about the S Matrix 
for M = 4 SYM at all loop level. As we will discuss in detail, the expression is also intimately 
connected to central ideas in algebraic geometry: Grassmannians, higher- dimensional residue 
theorems, intersection theory, and the Schubert calculus. This part of mathematics has not 
yet played a particularly central role in physics; it is both startling and exciting to find it 
sitting at the heart of scattering amplitudes in Yang-Mills theory. 

A given leading singularity is associated with a specific combination of residues or, equiv- 
alently, a choice of integration contour. In this paper we identify the contours corresponding 
to the leading singularities for all NMHV amplitudes and also the 8 pt N 2 MHV amplitude, 
leaving a complete exploration of the leading singularities for general amplitudes to future 
work. In the cases we consider, different representations for the amplitude correspond to 
different contours that can be smoothly deformed into each other: the remarkable identities 
guaranteeing the equality of all these representations, which at tree level also enforce the 
absence of unphysical poles, follow directly from the higher-dimensional analog of Cauchy's 
theorem known as the "global residue theorem" . The relations implied by the 1-loop "infrared 
equations" also follow directly from the global residue theorem. 

We believe that our conjecture represents a first direct look at the dual theory for the 
S Matrix that we have been seeking. The way in which such a simple formula manages to 
reproduce the incredibly complicated expressions appearing in scattering amplitudes is very 
striking-they are clearly being computed in a completely different way than local quantum 
field theory! Further explorations into its properties should lead to a more complete and 
physical formulation of the duality. With this goal in mind, we will study the properties of 
equation (Q with the aim of understanding what it is trying to tell us about M = 4 SYM. 

2 The Proposal 

Let us begin by defining some notation. We are interested in the single-trace n particle color 
stripped amplitudes in Af = 4 SYM, 

M n (X a , Xa,fj a ), a = l,---,n. (11) 
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Here we have chosen to label each particle by a \fj) Grassmann coherent state [4"4"t [9], where 
the 1 77 = 0) state is |+). The amplitude decomposes into the sum 

n 

M n = Y,M n , k (12) 

k=0 

where the M k have charge 4k under the U(1)r symmetry where fj have charge 1. As is 
well-known M k =o,i and M k = n ,n-i vanish; M n]k are precisely the N fc_2 MHV or, equivalently 
the N n_2 ~ fc MHV amplitudes. We will suppress the trivial dependence on the gauge coupling 
in everything that follows. 

Let us return to our object of study, this time giving it a name: 

(12-it)(23-(i t+ l))"..(nl-( t -l)) n^(g-W.) (13) 

where as before the W a are twistor variables obtained by fourier-transforming with respect 
to the A a , so W = {W\fj) = (/i, \\fj) and 

(mi • • • m k ) = e ai - ak C aimi ■ ■ ■ C aknik (14) 

Thinking of C aa as a k x n matrix, (mi ■ • • rrik) is the determinant of the k x k matrix made by 
only keeping the k columns mi, - ■ • ,m k . We'll use standard terminology and refer to these 
as the "minors" of C. 

We use the notation C n . k both because our claim is that this object computes the 'Xeading 
singularities" , as well to emphasize its relation to the "£ink" representation. Note that the 4k 
Grassmann 6 functions ensure that this expression has R-charge 4k as needed. The expression 
is manifestly superconformally invariant. It is manifestly cyclically invariant; in a moment 
we will also see that it is also parity invariant. However, at the moment the integral is rather 
formal and needs a proper definition. 

2.1 "Gauge Fixing" GL(fc) 

To begin, notice that the product of <5 4 ' 4 functions impose k linear relations on the n W a 's. 
These 6 functions are invariant under a GL(k) transformation 

C aa -> LiCfia (15) 

where L£ is any k x k matrix. 

A quick note before we proceed further: ultimately, we will be interested in complex 
kinematical variables, and our discussion will be valid for any space-time signature. However, 
for the moment we will keep things simple by defining our integral with real variables in (2, 2) 
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signature. In order not to clutter our expressions with explicit references to real and complex 
variables, we will write GL(fc) instead of GL(fc, R) or GL(k, C). Similarly, we will refer to 
the 4D Lorentz group as SL(2)x SL(2). 

It is very natural for the integral to be invariant under a full GL{k) symmetry, while the 
objects (mi • • • m k ) are invariant under the SL(k) subgroup. With a total of n minors in the 
denominator, the integrand weighted with the d kxn C measure has the full GL(k) symmetry. 
This completely fixes the form of the integrand up to our choice that the (m 1 ■ ■ -m fc ) consist 
of cross products of consecutive C 's of the form (j + ■ • ■ (j + k — 1) ). This seems to be the 
most natural choice, it is also quite likely that this form is dictated by dual superconformal 
invariance, though we won't pursue this further here. 

The GL(k) symmetry is like a gauge symmetry that generates a divergence in the d kxn C aa 
integral; to make sense of the integral we have to fix this gauge freedom. We do this by 
interpreting the C aa as a collection of n k— vectors. Using GL(k) we can set k of these 
vectors to any fixed set of k vectors we like. For instance we can set them to an orthonormal 
basis of the form (1, 0, • • • , 0), (0, 1, • • • , 0), • • • , (0, • • • , 0, 1). That is, we can choose a set 
bi, ■ ■ ■ , bk, and fix 

Cab fj = <W (16) 

Said another way, if C aa is a k x n matrix, then we are fixing k of the columns to an 
orthonormal basis of k vectors. This fixes k 2 of the C's, leaving us k x (n — k) free variables. 
We index the k particles belonging to the set bp with the letter /, and the remaining (n — k) 
particles with the letter i. The free C aa can then be written as cu- To be explicit, consider 
the case n — 7, k — 3. Here C starts as a 3 x 7 matrix, and we can gauge fix any three of 
the columns to an orthonormal basis. For instance one possibility is 



C 



( 1 Cu Ci 5 Ci 6 C17 

1 c 24 c 25 c 26 c 27 \ (17) 
\ 1 c 34 c 35 c 36 c 37 



while another is 




C = c 32 1 c 34 c 36 c 37 (18) 



Other gauge fixings of the GL(fc) may also be convenient, however in this paper we will 
canonically gauge-fix in this way. 

We can now define the "gauge fixed" k x (n — k) dimensional integral 



£ -< W "» - I (12.. .*)(23.. .(*+!)). ■(„!■■ .(fc-^ n^W + O.VV,) (19, 
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where the (mi • • • nif.) are computed from the C aa in one of its gauge fixed forms. 

Note that in this representation it is trivial to go to a basis where the particles indexed 
by / are in the Zj representation. Simply by fourier transforming, we find a form of the link 
representation [18] 



£ "*< W " Z ' ] = J (12-*)(23.. ■(*+!))■ .(nl-(t-l)) < 2 °> 

Thus, going from a basis in which k of the legs are in the Z representation to a different basis 
in which a different set of k legs are in the Z representation simply amounts to gauge fixing 
the C matrix in different ways. For instance the gauge fixing in equation (JTTj) naturally 
corresponds to the link representation of the amplitude with Z\, Z2, Z 3 , W4, W5, We, WV, 
while the one in equation ([TBI corresponds to Zi, W2, Z 3 , W4, Z 5 , W 6 , W7. 
With this form we can straightforwardly go back to momentum space: 

L. n . k {A,A, V) - J (i 2 ...jfc)(23...(A; + l)) ■■■(nl---(k-l)) [ ] 

Let us count the number of integration variables after the delta functions have been used 
to fix as many of them as possible. We have k x (n — k) link variables, and a total of 
2 (n — k) + 2 k — 4 independent bosonic 5 functions, where the —4 appears because four of 
these 6 functions ultimately become 4-momentum conservation. Therefore the number of 
remaining free integration variables is 

k(n - k) - (2n - 4) = (k - 2) x (n - 2 - A;) (22) 

Note that for k = 0, 1, n — 1, n the number of free integrations is negative. This reflects the 
fact that in these cases there are additional 5 functions in momentum space beyond the usual 
momentum-conserving one. For k — 0, the amplitude is just proportional to the product over 
all the S 2 (X)'s and so vanishes for generic momenta. For k = 1, the <5 2 's force all the A's 
to be proportional and this also vanishes for generic momentum-conserving momenta, as is 
familiar for the three-particle amplitude. Exactly the same argument holds for k = n, n — 1 
by reversing the roles of A, A. 

We can now state our conjecture precisely. As we have seen there is a (k — 2) x (n — k — 2) 
dimensional plane in en space, which are solutions to the equations 

A, - c n Xi = 0, A x + cnX, = (23) 

The solutions of the these linear equations for the en can be parametrized by variables r 7 , 
such that 

c/i(r 7 ) = c* H + ci/ i7 r 7 , with 7 = 1, • ■ • , (k - 2) x (n - k - 2) (24) 
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and where c* H is any particular solution to equation (123|) . Of course we can make any choice 
for the parametrization of the r 7 's that we like, for instance we can pick the r's to be any 
(k — 2) (n — k — 2) of the en, and solve for the rest of the c's in terms of them. 

Now, we can explicitly pull out the momentum conserving delta function from the 5 2 
factors as 

5 2 (X t - cnX^iXi + cn\ t ) = <5 4 (]>> a ) x J(A, A) x I d^^-^r^cn - Cji (r 7 )) (25) 

In this expression we cavalierly ignore the sign factors that arise in the Jacobian of real S 
functions; the sign would only appear as an overall factor depending on n, k and the external 
A, A in any case; we drop it since keeping it would introduce non-analytic factors not present 
in the amplitudes. In the appendix we present a more detailed treatment in which we discuss 
the relationship of our conjecture to the link representation of [18], and demonstrate that 
ignoring these sign factors gives the correct answer, justifying our loose treatment here. 
Our fully gauge-fixed and well-defined proposal is now 

£n-,k = L n . k x 5 4 C^2p a ) (26) 

a 

with 

r d (k-2)x(n-k-2) T 

Ln < k = J J [(12 ■ ■ • k) (23 ■ ■ ■ (Jfe + 1) ) ■ ■ • (nl ■ ■ • (k - 1) )] (r) H 6 ^ + (2?) 

At this point we can freely complexify the A, A and r variables. Since the integrand is 
holomorphic in the r's we must treat it as a contour integral in many complex variables. 

If we are interested in extracting the gluonic components of this supersymmetric object, 
it is particularly convenient to "gauge fix" so that the indices / which have been set to the 
orthonormal basis correspond to the negative helicity particles. Then, to obtain the gluon 
component we simply set fji —>■ to get positive helicity gluons and integrate J d 4 fjj for 
negative helicity gluons, and the Grassmann integral just gives 1. As an example, for the 7 
particle NMHV amplitude with helicities l + 2 + 3"4 + 5 - 6~7 + , we have n = 7, k = 3 and 

f d 2 r 
L ^ + + = 3 J [(123) (234) (345) (456) (567) (671) (712)] (r) (28) 

with the gauge- fixed C of the form 

/ c 3] 

C = C 5 l C 52 C54 1 C 57 I (29) 



C31 


C32 


1 


C34 








C37 


C51 


C52 





C54 


1 





C57 


C61 


C62 





C64 





1 


C67 
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2.2 Geometric Picture 



The momentum space form of our conjecture in equation fl2T|) has a very nice geometric inter- 
pretation. Indeed, we can motivate the conjecture by starting with an elementary observation 
about momentum conservation. 

Let us work in the n-dimensional "particle" space. The kinematical data is specified by 
giving X aa ,^aa- In any given Lorentz frame, we can think of X aa as simply labelling two 
n- vectors, A a= i, X a =2', similarly the A's correspond to two n- vectors Xa=i, Xa=2- Since Lorentz 
transformations act as SL(2)x SL(2) we should really think about the 2-plane A in the n- 
space, spanned by A a= i j 2, and similarly the 2-plane A. Note that the A, A planes intersect 
at what we can call the origin in n-space. Now in this n-dimensional setting, momentum 
conservation, J2 a ^aaXaa = is simply the statement that X a ■ A„ = 0, or put another way, 
that the A and A planes are orthogonal to each other. 

This is a quadratic constraint on the A, A. It is possible to "linearize" this constraint 
by introducing an additional auxiliary object. Consider a /c-plane C passing through the 
origin in the n-dimensional space. Associated to C is a natural (n — k) plane, C, which is 
simply the orthogonal complement of C. Given these objects, momentum conservation can 
be enforced by a pair of linear constraints: that C is orthogonal to A, and C is orthogonal to 
A. To spell out the obvious, if C is orthogonal to A, then C must contain the A plane, and 
since C is forced to be orthogonal to A, these two constraints together enforce momentum 
conservation. Note that these constraints are clearly impossible to satisfy for k = 0, l,n — 1 
or n, where either C or C is a point or a line. This gives a pretty geometrical explanation 
of why these amplitudes vanish. 



Let us now see how these geometric conditions are reflected in the equations. A k plane 
in n dimensions can be specified by a collection of k n-dimensional vectors, whose span give 





X i- P \ 
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the fc-plane. These k vectors can be grouped together into a k x n matrix C 



C 



( G x \ 

\c k J 



I C n 

C21 



C12 
C22 



C*2fc 



Ci,fc+i 

C2,fc+1 



C2,fc+2 



C*2.n-1 



C*l,n \ 
C*2,n 



k,2 



k,k+l 



c, 



k,k+2 



c, 



k,n—l 



c, 



(30) 



k,n J 



Naturally, this is a highly redundant description of a /c-plane since any k x k linear trans- 
formation on the k vectors leaves the fc-plane invariant. This is precisely the GL(k) "gauge 
symmetry" we encountered in the previous subsection. As before, we can "gauge-fix" the C 
matrix to put any k of the columns to some fixed set values, for instance 



C 



( 1 ... Ci jfe+ i C 1:k +2 
1 ... C 2 ,fc+l C 2 h+2 



... 1 Ckk+1 Ckk+2 



Cl,„-1 Ci 



C2,n-1 C2, n 



Ck,n— 1 Cfc,n J 



(31) 



and we see as before that /c-planes are specified by n x k — k 2 = k x (n — k) parameters. 

The space of fc-planes in an n-dimensional space is known as the Grassmannian G(k,n), 
and we have just seen that G(k, n) has dimension k x (n — k). There is a clear k (n — k) 
symmetry here, which reflects the natural association between the fc-plane and its orthogonal 
complement (n — A;)-plane. For the choice of "gauge fixing" made above, the in — k) plane 
C is given by an (n — k) x n matrix 



C 



( — Cl,fc+1 
— Cl,k+2 

\ ~ C l,n 



-C2,fc+1 
~C2,k+2 

~ c 2,n 



-Ck,k+l 
-Ck,k+2 



1 
1 



-Chn 







(32) 



Note that — signs in C are needed to ensure that C ■ C 
satisfied by the cu then take the simple form 

(7-A = 0, C-A = 



0. The linear equations 



(33) 



We can think of the representative C given above as defining a set of co-ordinates on the 
Grassmannian; in fact it is clear that this is only one chart covering an open set of all possible 
fc-planes, and that different choices of "gauge-fixing" the matrix C correspond to different 
charts covering different open sets. As we pointed out in the last subsection, in order to 
look at pure gluon amplitudes with a given helicity configuration, it is most convenient to 
gauge fix so that the columns corresponding to the negative helicity gluons are gauge-fixed 
to an orthonormal basis. This leads to a very pretty picture: the amplitude for different 
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helicity configurations correspond to integrating the same function, but over different charts 
on the Grassmannian! Keeping the supersymmetric Grassmann variables manifest instead 
ensures that we get precisely the same answer for the superamplitude no matter which chart 
is chosen. We will see how this works explicitly in some examples later. 

We have focused on momentum conservation in this discussion, but the superpartner of 
momentum conservation also follows explicitly from this geometrical picture. The amplitude 
is proportional to the Grassmann 5 function <5 4 (Cf/); since C must contain the A plane, we 
can pull out an overall factor of S 8 (J2 a ^aVa), which is precisely the superpartner of the usual 
momentum-conserving 5 function. 

In the last subsection we saw that after integrating over the bosonic 5 functions, we are 
left with (k — 2) x (n — k — 2) free variables. Not coincidentally, this is the dimension of 
G(k — 2, n — 4), as is completely obvious geometrically. Recall that we are constraining C 
to be orthogonal to A and contain A. Since C must contain A, of the k vectors needed to 
define C, it is natural to choose two of them to span the plane A, and choose the remaining 
(k — 2) to be orthogonal to both A and A. Thus the /c-planes satisfying our constraints are 
naturally in one-to-one correspondence with a (k — 2) plane in the (n — 4) dimensional space 
orthogonal to A and A. We can think of the C(r 7 ) matrix satisfying the constraints as a 
mapping from G{k — 2, n — 4) — > G(k, n), which we are integrating with a natural measure. 

This finally brings us to the interpretation of the minors (i, i + 1 ••• ,i + k — l) appearing 
in the integration measure of our formula. Thus far no particular basis in the n-dimensional 
space has played a privileged role in our discussion, however, the external states and color 
ordering of do make the orthonormal basis set e 1 = (1, 0, • • • , 0), • • • , e n = (0, 0, • • • , 1) 
special. The k n-vectors defining the /c-plane can clearly be projected into the A;— dimensional 
subspace e mi , • • • , e mk and will generally fill some volume in this space; the minor (mi • • • rrik) 
is that volume. By itself this volume is not a particularly natural geometrical object, since 
it is not invariant under the full GL(k) symmetry but only under SL(/c). However, as we 
commented in the last subsection, due to the transformation of the measure, the whole 
integral is nicely GL(k) invariant. Indeed, this motivates the interpretation of the set of 
all (™) minors {m 1 ---m k ) as coordinates in a pW -1 projective space. This embedding of 
G(k, n) into pW^ 1 is known as the Pliicker embedding, and the minors are referred to as the 
Pliicker coordinates. 

2.3 Manifest Cyclic and Parity Symmetries 

Before moving on to some simple examples, let us quickly show that the cyclic symmetry 
and parity are manifest in our proposal. By "are manifest" , we really mean that they leave 
the integrand invariant; of course the choice of contour can break these symmetries. 
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First for the cyclic symmetry: this is trivially present in the original form of the proposal. 
While the gauge fixing breaks the manifest cyclic invariance, the underlying gauge symmetry 
ensures that all the different gauge fixings yield the same cyclically invariant result. We will 
see this explicitly in a number of examples. 

While the cyclic symmetry is obvious before gauge-fixing, parity is obvious after gauge 
fixing. Our integral is 

£ - = I (12... fe )(23..^ C .(nl... fe -l) ^ ' ^ ' ' '<> < 34 > 

Fourier transforming with respect to fj and swapping k «-> (n — k) amounts to also swapping 
C <-> C. It is also easy to see that in any gauge fixing, the minor (12 • • • k) of C is the same 
as the minor (k + 1 ■ ■ - n) of C, and thus the measure and integrand are parity invariant. To 
see this in a concrete example, consider the 6 particle MHV amplitude. Gauge fixing the 
first two columns of C, we find that C and C are 

C=(\ ? Cl3 ° U Cl5 Cl6 ) (35) 



C = 



1 c 23 c 24 c 25 c 26 

/ -ci3 -c 23 1 \ 

cm -c 24 1 

cis -c 25 1 

V -C16 -C26 1 / 



(36) 



and we see immediately that (45) = (6123) = Ci 4 c 25 — Ci 5 c 24 , and so on. The general statement 
is easily proven by induction. 

2.4 The Polynomial Degree of the Minors 

Let us make one final general observation, and determine the order of the polynomials in the r 7 
that will appear in each determinant (mi, • • • , m^). In general the solutions cnij) = c* H + Tn 
where th = d/i 7 r 7 are linear in the r 7 and satisfy 

T/;Aj = 0, TnXi = (37) 

The first equation tells us that, thinking of set of Ti /c-vectors, all these vectors are 

orthogonal to the 2-plane A. Therefore all the r lie on a (k — 2) dimensional subspace, and 
so the cross product of any (k — 1) of them must vanish. Hence, in the k cross-products 
appearing in (mi • • -m^), the th variables can appear at most (k — 2) times, and each factor 
in the denominator is therefore a polynomial of degree at most (k — 2) in the free integra- 
tion variables. Note that we could have made exactly the same argument using the parity 
conjugate form of the amplitude which would tell us that the polynomial is of degree at 
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most (n — k) — 2. Thus in general each factor in the denominator is a polynomial of degree 
min[(£; - 2), (n - 2 - k)\ in (A; - 2) x (n - 2 - k) variables. The NMHV amplitudes are 
particularly simple: there are (n — 5) integration variables and each term in the denominator 
is linear in them. 

3 First Examples 

3.1 MHV Amplitudes: k = 2 

In this case we set k = 2. From our general formula we can see that there is no integration to 
be done and therefore L n . 2 is straightforward to evaluate. Before evaluating L n . 2 in detail, let 
us use this example to illustrate some of the geometrical properties we discussed in generality 
before. The geometrical ideas provide an intuition for motivating our formula. 

3.1.1 Direct Geometrical Evaluation 

In section 2.2, momentum conservation was expressed as a fully geometrical condition. Two 
fixed 2-planes, the A and A plane were introduced. Momentum conservation is simply the 
statement that these two 2-planes are orthogonal in C n . 

For k = 2, we are interested in the space of all 2-planes that contain the A-plane and 
are orthogonal to the A-plane. Clearly there is only one such plane and it is must coincide 
with the A-plane itself. This means that the 2 n-vectors giving the rows of the 2 x n matrix 
representation of C must be linear combinations of (A|, \\ , . . . , A") and (A2, A|, . . . , A??). Since 
we are interested in the determinants of the minors, we are free to choose them to be exactly 
equal to the two A n-vectors, 

n _ ( C11 C12 C13 • • • Cl,n-1 Ci,n \ _ ( Aj A^ A^ ... A™ 1 A™ A , , 
\ r r r n r I ~ \ \ l \ 2 \ 3 x™^ 1 \ n I ' ^ ^ 

V <^21 <-^22 <-^23 • • • ^2,n-l < - / 2,n / \ ^2 A 2 A 2 • • • A 2 A 2 / 

This fixes the GL(2) "gauge symmetry". With this identification it is clear that the determi- 
nant of a 2 x 2 minor made from columns i and j gives (ij) = (ij). Therefore, the product 
(12) (23) ... (n — In) (n, 1) gives rise to the usual Parke- Taylor formula for the denominator 
of an MHV amplitude. 

The numerator is fixed by imposing the supersymmetric version of the geometric con- 
dition: C must be a 2-plane which is orthogonal to a 0|4 plane in C n '" spanned by the 4 
n-vectors (fj(, fj^, ■ ■ ■ , fj^), which is imposed by the Grassmann 5 functions 

n n n 

a=l a=l a=l 
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Putting everything together we find 

Ln ' 2 - (12}(23}...(nl)- (40) 

In this derivation we have made an identification which only works for k = 2; we have 
embedded the SL(2) action of the Lorentz group inside the GL(2) acting on the space of 
2 x n matrices. 

3.1.2 Evaluation Using Canonical Gauge Fixing 

Let us now show how the same formula can be obtained by the canonical gauge fixings of 
the GL(2). 

Put k = 2 in equation (|2ip . and gauge fix so that the index / runs over the particles x, y. 
Therefore, on the support of S 2 (Xi — cnXi), we can solve 

Xj = cnXj =► c u = (41) 

(xy) 

Note that we can more generally write 

X a = C aa X a =>- A a A b = (xy)e^C aa C pb = (xy)(ab) (42) 
where the a, (3 indices also range over x, y. It is also trivial to see that 

5 2 (X - cuX^CXj + CI X) = J5\Y,p^ with J = ^y) A ' n ( 43 ) 

a 

Therefore 



1 (xyY 



(12)(23)---(nl) (12)(23)---(nl) 
The Grassmann delta functions can also be easily simplified starting with 



(44) 



) (45) 



and using (141]) we obtain 



1 n 

<^ 8 E A <^) ( 46 ) 



a=l 

Combining (jUJ) with fT4oT) we get the desired result. 

Note that already with this simple example we see explicitly something that we claimed 
on general grounds: the final form of the amplitude is independent of how we "gauge-fixed" 
the GL(2) symmetry. 
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3.2 6 Particle NMHV Amplitude 

Next, let us consider NMHV amplitudes with k — 3. The number of integration variables is 
(3 — 2) x (n — 3 — 2) = {n — 5). Obviously the simplest case is with n — 6, which involves 
only one integration variable which we call r. Since this is such a simple case, we will work 
through the computation of the residues of interest rather explicitly here. 
The relevant integral for the 6 particle NMHV amplitude is 

/d 
[(123)(234)(345)(456)(561)(612)](r) ^ ^ + (4?) 

To study the alternating helicity amplitude, we choose a convenient gauge fixing 

(c 2 i 1 c 23 c 25 \ 
C4i C43 1 C45 J . (48) 
C61 c 63 c 65 1 / 

In terms of the en appearing in this matrix, the minors are given by 

(123) = c 41 c 63 - c 6 ic 43 = c 25 , (345) = c 41 , (561) = c 63 , 

(234) = c 36 , (456) = c 25 , (612) = c 41 . (49) 



so that 

1 1 



(123) (234) (345) (456) (561) (612) c 25 c 41 c 63 c 25 c 41 c 63 

and 

dr 



^6,3 



(50) 



(51) 



[c 2 5C 63 C4iC 2 5C 63 C 41 ](r) 

Here we have introduced the notation en to denote a pole which maps to en, and vice- versa, 
under a parity transformation. 

To identify cu(t) explicitly, we are looking for the 1-dimensional space of solutions to the 
equations 

\ - cnXi = 0, A/ + c n \i = (52) 

Since any three two-dimensional vectors are linearly dependent, they satisfy the "Schouten 
identity" which we write in the form 

e ijk \\jk} = 0, euKXjiJK) = (53) 

Hence, given any particular solution c* H , we can find another solution 

cr(t) = c* H + e ijk e IJK [jk]{JK)T (54) 
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With this choice it is easy to verify that the Jacobian J is 1: 

n ^ - n 52 & + c «^) = 1 x a ^«) / dr69 - ( 5s ) 

i I a J 

We can also always choose the origin for r so that, for example, 025(0) = 0; this allows us to 
solve for the c* H . For instance we can use the A 5 and A 2 equations to solve for c 65 , c 45 , c 2 i, c 23 , 

e - (54) c* - (56) c* - [23] c* - [21] f56l 
c 65 - (64) , c 45 - (4g) , c 21 - [13] , c 23 - [31] (5b ) 

and then use, say, the A 4 , A 6 equations to solve for the rest of the c's 

(6|(5 + 4)|3] _ (6|(5 + 4)|1] (4|(5 + 6)|3] (4|(5 + 6)|1] 

41 (46) [13] ' 43 (46) [31] ' 61 (64) [31] ' 63 (64) [13] 1 J 

In the above computation we have used a tiny bit of foresight to parametrize the one- 
dimensional space of solutions, but we would have arrived at precisely the same expression 
by brute force, for instance by picking one of the c's, say c 2 5, to be special and solving for 
all the rest of the c's in terms of this one. That would be a different parametrization of the 
1-dimensional space of solutions with a different J, but of course precisely the same final 
answer. 

As we saw on general grounds, each of the minors is linear in r, and so the integrand 
has 6 poles, associated with the points where each minor vanishes. Let us introduce some 
compact notation to denote the relevant residues. A given minor (i, % + 1, % + 2) is specified 
by its starting point %. We will therefore refer to the residue at the pole corresponding to the 
vanishing of the minor (i, i + i + 2) as {i}. 

A short computation then yields 

{1} = " [23][34](56>(61>(p 5 +1 + p!)=<l|6 + 5|4](5|6 + 1|2] ' {3} = ^ {1} ' {5} = 9 ' {1} (58) 
and 

{4} = [12] [23] (45) (56) (p 4 + p 5 + J 2 (6|5 + 4|3](4|5 + 6|1] ' W = ^ 4 >' ^ = »^ 4 > ^ 

The reader will recognize the residues {2}, {4}, {6} as the three terms in the BCFW form 
of the 6 particle tree amplitude, and the {1}, {3}, {5} residues as the negative of the three 
terms in the P(BCFW) form of the same amplitude: 

A4cfV _ = {2} + {4} + {6}, M p + (B + F +y = -{1} - {3} - {5} (60) 

We can now identify the tree amplitude with a particular choice of contour in our integral; 
depending on what contour is chosen we can get different forms of the amplitude. Of course 
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since we have the full fj dependence, we can obtain any helicity amplitude we please, but for 
simplicity let us start with the alternating helicity amplitude which is naturally associated 
with our form of the gauge-fixing. 

Let T be a contour that encircles the poles where (234), (456), (612) vanish, and F be 
a contour that encircles the poles (123), (345), (561). Note that under a cyclic shift, the 
poles contained in T map to the ones contained in T and vice-versa. Most naively, then, 
to extract a cyclically invariant object from the residues we should take the contour T + T; 
in fact there is a minus sign in the parity mapping from T to T and the correct, manifestly 
cyclically invariant contour is T — T. This gives the manifestly parity symmetric form [BCFW 
+ P(BCFW)] of the amplitude. However, since the integrand vanishes as -% as r — > oo, up 
to a factor of 2 we can also write the amplitude using only the T or the — T contour. These 
give us the BCFW and P(BCFW) forms of the amplitude individually. Applying Cauchy's 
theorem to the contour T + T and enclosing all the poles with the same orientation gives us 
the remarkable 6 term identity which guarantees the equality of the BCFW and P(BCFW) 
forms of the tree amplitude. 




As we remarked in the introduction, the 6 term identity not only guarantees the cyclicity 
and parity of the tree amplitude, but also enforces the absence of unphysical poles. This is 
also easy to see from our contour integral. As we change the external kinematics, the position 
of the 6 poles in the r plane move; we should only expect singularities when poles collide. 
The "unphysical" poles correspond to a collision between poles contained in e.g. T, or ones 
in T. For instance, looking at the residue where (456) = C25 = 0, when (6|5 + 4|3] — > 0, the 
minor (612) = C41 — > 0. However clearly there is no actual singularity here, since we can 
always deform the contour to — T, none of whose encircled poles are colliding. The physical 
singularities, on the other hand, involve the collision of a pair of c and c poles; since the tree 
amplitude contour separates the c poles from the c poles, the contour is necessarily pinched 
between them and a singularity arises. 
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Let us look at the physical singularities in more detail. Given that there are 3 c and 3 
c poles, there are 9 ways of colliding them in pairs. On physical grounds we expect poles of 
the form si 2 , S23, S34, S45, S56, sqi, £123, £234, £345- Note that collinear poles can be factored 
sa + \ = + [ii + l] while ta+i i+2 cannot. This means that t-poles must be parity invariant. 
Indeed, note that a £123 pole is reached when (123) = c 2 5 and (456) = c 2 5 collide. The same 
happens for £234 and £345. Working with complexified momenta means that we can consider 
the poles (ij) and [ij] as independent. This gives 12 poles but there are only 6 pairs of 
non-parity related poles that can collide pinching the contour. This means that complex 
collinear limits must also pair up. It is easy to see that this is the case. For example, poles 
c 2 5 and c 4 i collide when [12] — > or when (45) — > 0. 

Let us quickly see how to obtain the other helicity configurations. For instance, consider 
the "split-helicity" configuration M h++ . Performing the appropriate fj integrations sim- 
ply multiplies the alternating helicity integrand by a factor C25(t) 4 . This has the effect of 
removing the pole at C25 = 0, though the integrand still vanishes at infinity like Thus the 
BCFW form of this amplitude has only two terms, while the P(BCFW) still has three. 

It is also instructive to see how all of these results can be recovered from another gauge 
fixing. For instance, let us consider the gauge fixing where the first three columns are set to 
the orthonormal basis 

/ 1 c u ci 5 cie \ 
C = 1 c 24 c 25 c 26 (61) 
\ 1 c 34 c 35 c 36 / 

This is the gauge- fixing convenient for the split-helicity configuration; we see that (123) = 1, 
(234) = C41, (345) = C14C25 — C24C15, and so on. From this form it is immediately apparent 
that, in computing the split helicity amplitude, we only have 5 poles in the complex r plane, 
since by the gauge-fixing the minor (123) is identically set to 1. But how are we to see all 
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six terms needed for the alternating helicity amplitude? The answer is obvious. If we want 
to use this gauge-fixing to compute the alternating helicity amplitude, the appropriate fj 
integrals once again multiply the integrand by C2s(r) 4 . This grows as r 4 for large r, while 
the denominator falls as r -5 , so there is a pole at infinity; this precisely gives the "missing" 
residue, which we have checked correctly completes the amplitude. 



3.2.1 All Loop Leading Singularities? 

There is something a little peculiar about our discussion of the 6 particle NMHV amplitude. 
The residues of our object appear to contain more information than what is needed for 
the tree amplitude, which are a particular combination of the residues. It is natural to wonder 
whether the individual residues have any meaning. As we reviewed in our discussion of one- 
loop leading singularities, BCFW and P(BCFW) terms are the 1-loop leading singularities 
associated with scalar boxes with two adjacent massless legs. But all the scalar boxes with 
6 external legs are of this form. Therefore, all the leading singularities of the 6 particle 
NMHV amplitude are determined by the residues. The box coefficients are given by residues 
as shown below 
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where all other boxes are cyclically related to thes^l There are two terms since each box 
receives a contribution from a leading singularity and its parity conjugate. 

With this in mind, we can revisit the case of the MHV amplitude, where L n;2 involves 
no integrations and directly produces the tree amplitude. The one-loop MHV amplitudes in 



2 The fact that the coefficient of these four scalar boxes is the same guarantees the IR equation from the 
double- logarithm of £123- 
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M = 4 were computed in [45] and are given by 

^C;MHV = K;MHV IP,Q>R,S ( 62 ) 

{P,Q,-R,S}e{lm,2m-e} 

where the sum is over all one-loop 1-mass and 2-mass-easy box integrals normalized to have 
leading singularity equal to one. Thus, the only non-vanishing leading singularities are pro- 
portional to the MHV tree amplitude, and one can therefore also give a 1-loop interpretation 
to L n . k for MHV amplitudes! 

So far the evidence that L n -k is actually computing 1-loop leading singularities is rather 
circumstantial; we will see more dramatic and direct evidence starting with the 7 particle 
NMHV amplitude. But first, let us take a peek at our ultimate claim, which is that L n .^ 
computes leading singularities at all loop orders. In particular, all leading singularities of 
two-loop amplitudes have been computed for MHV amplitudes up to 6 particles and have 
been found to either vanish or be proportional to the tree MHV amplitude. The reader might 
wonder how a two-loop computation can possibly give rise to such a simple object as a tree 
amplitude. In order to illustrate this we consider the four-particle 2-loop amplitude and show 
that it provides the prototype of what we called "composite" leading singularities in section 
1.1. 

Following the discussion and the notation in section 1.1., at two-loop level and four 
particles, the closest we can get to 4 x 2 = 8 propagators is to collect all Feynman diagrams 
with seven propagators with a "double-box" topology. After choosing the contour where all 
seven propagators are on-shell, the rational function R(£i,£ 2 ) in 

/ t~t7 pi R \^*V ( 63 ) 
■> 1 li=l i 

factors as the product of seven tree amplitudes shown in the figure. The piece on the right 
hand side is a one-loop leading singularity which is equal to a 4-point tree-amplitude. This 
four-point tree amplitude has two factorization channels. The remaining loop integral can 
be taken as a contour integral on a T 4 defined so that the t-channel of the four-particle tree 
amplitude is manifest. We are then left with a 1-loop leading singularity which evaluates to 
a four-particle tree amplitude once again! This is precisely the procedure used to obtain the 
5 and 6-point results. 




26 



There is also a computation in progress for the 6 particle NMHV amplitude |43j . Quite 
remarkably, for the 6 particle NMHV amplitude, the very same rational functions we have 
already encountered at 1-loop determine the 2-loop leading singularities as well! Thus, the 
available data at higher-loop orders is non-trivially consistent with the idea that L n .^ is indeed 
computing leading singularities at all loop orders. We will provide much more compelling 
evidence that two-loop leading singularities are present in our discussion of the 8 particle 
NMHV and N 2 MHV amplitudes. However, let us first consider the striking evidence for 
1-loop physics in the 7 particle NMHV amplitude. 

3.3 7 Particle NMHV Amplitude and Direct One Loop Evidence 

Consider k = 3 and n = 7. The number of integration variables is 2. From the six particle 
example we learned how different gauge fixings lead to equivalent forms for the amplitudes. 
We find it most useful to consider "gauge" choices in which there are no contributions from 
poles at infinity. 

In this discussion we concentrate our attention on the phenomena that do not appear for 6 
particles. The most important is the definition of residues in more than one complex variables 
and the generalization of the residue theorem which we will find leads to the equations 
obtained from the IR behavior of one-loop amplitudes. 

To be concrete let us look at the (l~2 + 3~4 + 5~6 + 7 + ) amplitude. Then we have 

r d 2 r 

Lr ' 3 ~ J J [(123) (234) (345) (456) (567) (671) (712)] (r) (64) 
where the matrix C is of the form 

(1 C12 C U Ci 6 Ci7 \ 
c 32 1 c 34 c 36 c 37 (65) 
c 52 c 54 1 c 56 c 57 / 

As we have mentioned in this case each of the factors (i, i + 1, % + 2) in the denominator is 
linear in the two variables (r 1; r 2 ). 

3.3.1 Multivariable Residues 

We want to interpret the integral as a contour integral. In order to learn how to compute 
residues in this situation, let us consider a simple example. Consider the following function 
of two complex variables Z\ and Z2 

f{ZU Z2) = (a^ + bzi + cKet+fzi + g) (66) 
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where h(zi, z 2 ) is any function which is regular where (azi+bz 2 +c) = and (ezi+fz 2 +g) = 0. 
Then we want to define a contour integral of the form 

j — f h(z 1 , z 2 )dzidz 2 ^ 

J (az 1 + bz 2 + c)(ez 1 + fz 2 + g) 

which can be called a residue of /. The most natural way to evaluate an integral of this form 
is to perform a change of variables U\ = az\ + bz 2 + c and u 2 = ez\ + fz 2 + g. Therefore, we 
have 

dui du 2 h(zi(u), z 2 {u)) 



I 



Ul U2 det(f^) 



(68) 



Now it is very natural to define the residue at (1x1,1x2) = (0,0) as the integral over = e 
and \u 2 \ = e for some small real number e in complete analogy with the one-dimensional case. 
Therefore the residue is given by 

Res[/](,r,, 2 *) = • (69) 

det ( ) 

Here z\ and z\ are the solutions to az\ + bz 2 + c = and ez\ + fz 2 + g — 0. 

Note that due to the antisymmetry of the determinant in (ui,u 2 ), the residue defined 
in this way depends not only on {z*,zl£) but also on the order in which the two vanishing 
factors Mi, u 2 are written. This is unlike what we are used to with a single complex variable; 
the reason is that in the familiar case, the contour defining the residue actually encloses the 
pole. In the case at hand, the contour of integration used in the definition of the residue is 
the product of two circles, S 1 x S 1 ^ T 2 in C 2 , which does not enclose the point (zl, z 2 ). In 
the math literature this T 2 is called the distinguished boundary to emphasize this fact. 

Using this result, let us consider a function of the form 

Pi{z u z 2 )p 2 {zi, z 2 ) ■ ■ •Pm{zi,z 2 ) 

For M > 2 and the Pi(zi,z 2 ) are linear. Then we have (^) residues, determined by putting 
any choice of two of the factors Pi^p^ to 0. The residue is given by 

Resf/K^*) = giZlZ; \^. • (71) 



rW)P*(^ 2 *) det (^>) 



We can trivially generalize this discussion to functions of iV complex variables of the form 

/ (zi,--- ,z N = — r 72 

Pi(z u ■■■ ,z N )-- -Pm\z u ■■■ ,z N ) 
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where the pi (z) are linear, and M > N. We can choose N of denominator factors, p% x ,- ■ ■ , Pi N , 
and solve the linear equations Pi x {z{, ■ ■ ■ , z* N ) = ■ ■ ■ = p iN [z\, ■ ■ ■ , z* N ) = to determine the 
point z*, - - • ,z* N . The residue is then defined to be 

Res[f)(zl ■ ■ ■ , z* N ) = 9{Zl ['' lZ * N ] a{Pi —^ (73) 

Ui^n.-^Pi^h ■ ■ • , Z*N) det ( d(zl,-,z£) ) . , 

Note again that, since the determinant is antisymmetric in (ii, ■ ■ ■ , ijv), the residue also has 
this antisymmetry property. So, the residue is not only associated with the point z*, ■ ■ ■ ,z* N) 
but the sign of the residue is determined by the order in which the factors Pi x , ■ • • , pi N are 
taken. 

3.3.2 One Loop Leading Singularities 



In equation (|64p . one has seven linear factors in the denominator. This means that we have 
(2) =21 different residues. We choose to denote each residue by the pair of linear factors 
which vanish; we denote the residue obtained by putting (i, i + i + 2), (J, j + 1, j + 2) — ► 
as {i, j}. Note that as mentioned above this residue is naturally antisymmetric, so that 

{ij} = -0,0- 

The computation of the residues is completely straightforward, and only involves solving 
linear equations. Just to illustrate this first non-trivial example with more than one complex 
variable we give some details in the appendix, but for the rest of this paper we leave the 
computation of the residues in the capable hands of Mathematica. Even though the residues 
only involve solving simple linear equations, the expressions quickly become complicated. To 
give two examples, we have 

r, = (13) 4 [7|(4 + 6)|5) 4 

1 ' 1 <12> <23> <45> <56> [7| (5 + 6) |4> [7| (1 + 2) |3> <6| (4 + 5) (2 + 3) 1 1) 

and 

([7|(2 + 4)|3>(54> + [76]<65><34» 4 



{2,4} 



(23) (34) (45) (56) [71] [1 1 (2 + 3) |4) [7| (5 + 6) |4) (4(5 + 6) (7 + 1) |2) <4| (2 + 3) (7 + 1) |6) 



We will spare the reader the sight of all 21 residues. 

The connection of these 21 residues to the corresponding one-loop amplitude is the fol- 
lowing. In [16], Bern et.al. computed the full seven-particle amplitude. In this case one has 
1-mass, 2-mass-easy, 2-mass-hard and 3-mass boxes. In [16] it was found that all the coeffi- 
cients could be written as linear combinations of basic objects. The amplitude ( — I 1 h+) 

has a flip symmetry which is obtained by applying i — > mod(6 — i) + 1 to the particle labels. 
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In [16], a list of 12 basic objects was given. Out of these, 3 are flip invariant. This means 
that the total number including the flip images is given by 3 + 2x9 = 21. Quite remarkably, 
our 21 residues map one to one to the 21 basic objects in [46] . 
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We can also specifically identify the residues that give the tree amplitude. Once again 
there is a BCFW and P(BCFW) form of the amplitude, and we can identify them amongst 
the residues as 



^bcfw _++ = {1,3} + {1,5} + {1,7} + {3, 4} + {3, 6} + {5, 6}, (74) 
M p"(BCF W r = -{2, 3} -{2, 5} -{2, 7} -{4, 5} -{4, 7} -{6, 7}. (75) 

3.3.3 Relations Between Coefficients 

The precise agreement between our 21 residues and the 21 rational functions first defined 
by Bern et.al. in the 1-loop amplitude is a quite remarkable confirmation that our object is 
computing 1-loop amplitudes. This agreement is even more surprising given the fact that, as 
noted in [16], these basic objects are not linearly independent. In [16] four relations among 
the 21 basic objects were given. One of them is flip symmetric. This means that counting 
the flip images one gets 1 + 2x3 = 7 relations. These relations come from the IR equations 
for the 1-loop amplitude, which as reviewed in the appendix, relate the double-logarithmic 
divergences of the 1-loop amplitude to the tree amplitude. These equations constrain the 
coefficients B of the scalar boxes; partitioning the external states into 4 sets (31,52,53,34), 
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the coefficient of the box with the sums of the momenta in S{ flowing into its corners is 
denoted by £( S1 )( S2 )( S3 )( S4 )- 

Let us look at the equation associated with the double- logarithm of £123, 

= — -B(l)(2)(3)(4567) + #(7) (1) (23) (456) + #(3) (4) (567) (12) ~ ^-8(4)(5)(67)(123) ~ | -8(6) (7) (123) (45) (76) 
--8(4) (56) (7) (123) + 5 (1) (23) (4) (567) + #(3) (456) (7) (12) + | -8(3) (45) (67) (12) + | -8(1) (23) (45) (67) • 

Here e.g. -B(i)(2)(3)(4567) is the coefficient of a box with momenta flowing through the corners 
given by the sum of momenta appearing in the parentheses. Using the explicit map we have 
just found between leading singularities and residues we find 

B(7) (1) (23) (456) = {1, 4} + {4, 5}, -B(3)(4)(567)(12) = {5, 1} + {4, 5}, -B(4)(5)(67)(123) = {5, 1} + {1, 2} 
5(6) (7) (123) (45) = {1, 4} + {7, 1}, #(4) ( 56 ) ( 7 ) (123) = { 4 , 5}, #(1) (23) (4) (567) = {1, 2}, 

8(3) (456) (7) (12) = {7, 1}, 5(3) (45) (67) (12) = {3, 1}, 8(i)( 2 3)(45)(67) = {1,6} 

and 

8(i) ( 2)(3)(4567) ={1,4} + {4, 5} + {1,6} + {1,2}. (77) 
Using this in equation (176|) one finds a surprising simple result 

= {1,2} + {1,3} + {1,4} + {1,5} + {1,6} + {1,7}. (78) 
Clearly, by cyclic symmetry we have seven relations of the form 

7 

= ]T{j,z}, j €{!,..., 7} (79) 

i=l 

Note that {i,i} = 0. 

In our particular problem, we have a flip symmetry. The relation coming from j = 2 is flip 
symmetric while the other six split into pairs related by a flip transformation. Very nicely, 
in |46j . three of the four relations given are six-term relations of the form (!79l) while the last 
one is a 12 term identity that can be shown to follow from the other six-term identities and 
one more of our six-term identities. 

Expressing the residues as explicit rational functions of the kinematical invariants, these 
identities appear to be miraculous statements. Another amazing identity, completely anal- 
ogous to the six-term identity we encountered with 6 particles, is the 12 term identity 

^BCFW = ^P(BCFW)- 

3.3.4 Relations Arise Prom Generalized Residue Theorems 

The form of the IR equations (!79|) begs for a residue theorem derivation. Note that we have 
seven relations. In the theory of residues of a single complex variables one expects a single 
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relation. This is the first difference between single and multidimensional residues. We will 
touch on these points in more detail and generality in the next section, but for the case at 
hand we can make do with repeated application of the usual Cauchy theorem. Our integral 
is of the form 

n , dTld I 2 r (80) 

lL( a i r l + b i T 2 + Ci) 

in our example the product is taken for i — 1, ■ • • ,7 but we'll imagine it for any i — 1, • • • , M 
for M > 2. Let us choose the j th factor in the denominator and treat r± as a fixed variable. 
Therefore we can use the integral in r 2 to compute the residue at r 2 = (— Oj-ti — Cj)/bj. The 
remaining n integral is then of the form 



I 



r 



IL/j&Ti + Ci) 

Now the usual one dimensional residue theorem tells us that the sum over the M — l residues 
is zero. Recalling that each residue came from first setting + 1, j + 2) = we find 

= £O,i}. (82) 

i 

As we saw in the previous subsection, these equations are nothing but the IR equations for 
the 1-loop amplitude! Thus, for the 7 particle amplitude, the IR equations which imprint 
both locality and unitarity in the 1-loop amplitude, are a direct consequence of the residue 
theorem for two complex variables. 

There is a simple corollary of equation (152"]) that is sometimes useful. Let if be a subset 
of {1, ... , M} and H its complement. Then the following set of equalities hold 

= LEW} = ££{m} = EEfttf ( 83 ) 

The last equality holds due to the anti-symmetry property. Consider the identity obtained 
by choosing H = {1, 2, 4, 6} and H = {3, 5, 7} for our seven-particle amplitude. The identity 
is 

= {1, 3} + {1, 5} + {1, 7} + {2, 3} + {2, 5} + {2, 7} + {4, 3} + {4, 5} + {4, 7} 

+{6, 3} + {6, 5} + {6, 7}. (84) 

It is easy to see that this 12-term identity is precisely the one obtained from equating the 
expressions for M B cfw an d M P ( B cfw) i n equation (1741) . In this example we appeared to 
need some foresight to make the correct choice for H to produce the needed identity for 
tree amplitudes. In fact there is a more natural and systematic way of understanding these 
identities, which we defer to a general discussion of all NMHV amplitudes in section HI 
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3.3.5 Prediction for Higher Loops 

We have seen that every one of the residues associated with the 7 particle NMHV amplitude 
can be identified with 1-loop leading singularities. In complete analogy with what we saw for 
the MHV and 6 particle NMHV amplitudes, we then predict that the leading singularities 
for the 7 particle amplitude at two loops and beyond should all be determined by the 21 
objects we have already identified. In the next subsection, we find more direct evidence for 
the presence of new objects associated with 2-loop leading singularities starting with the 8 
particle NMHV amplitude. 



3.4 8 Particle NMHV and Direct Two Loop Evidence 

The expressions for the box coefficients of general 1-loop NMHV amplitudes have been de- 
termined in [H], and we will discuss the identification of these with residues in the next 
section. But in this subsection we look at the 8 particle NMHV amplitude specifically, be- 
cause for the first time in this case not all residues are accounted for amongst the 1-loop 
leading singularities, and we can identify new 2-loop leading singularities. 

The mapping between box coefficients and residues is given below. We use a notation for 
residues which is the obvious generalization of what we introduced for the 7 particle case: the 
residue {i, j, k} is associated with putting the minors (i, i + 1, i + 2), (j, j + 1, j + 2), (k, k + 
1, k + 2) to zero, and is antisymmetric in exchanging any of the indices. For the 3-mass boxes 
and the 2- mass-hard boxes with two adjacent massless legs, we have 
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For the 2-mass-easy boxes with massless legs at opposite corners, and the 1-mass boxes, we 
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have 





/ 








/ 


\ 



7 



1 



H 3 



$1,2,3,] + $5,6/7] $1,2,h] + $1,2,6] + $1,2,8] 



7 2 



^1 / 
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3 a 

$3,f ,5] + $1,6 ,7] + $3,6,7] + $5,6,7] + $1, V] + $3,4,7] + $1,1,7] 

All other boxes are related to these by the cyclic symmetry. 

We defer a discussion of the tree amplitudes to the next section, where we discuss these 
for all NMHV amplitudes in generality. 



3.4.1 IR Equations and Residue Theorems 

Let us take a very quick look at some of the IR equations and their origin in residue theorems. 
Consider the IR equation associated with a double-logarithmic dependence on £1234. This is 
a particularly simple case to look at, since the 1-mass box, which has the longest expression 
in terms of residues, does not produce this singularity and therefore does not participate. 
This IR equations involves 14 box coefficients: 



2 5 (7)(8)(1234)(56) 



2 5 (1)(2)(34)(5678) 



2" B (3)(4)(5678)(12) ~ #(8)(1234)(5)(67) 



--S(l)(234)(5)(678) ~ #(1) (23) (4) (5678) + #(8) (123) (4) (567) + ^"^(l) (234) (56) (78) 



-B, 



(5) (678) (12) (34) 



-B, 



(8) (12) (34) (567) 



+ -B, 



(4) (56) (78) (123) 



+ B 



(4) (5) (678) (123) 



+B 



(8) (1) (234) (567) 



2 -6(5) (6) (78) (1234) = 



J5) 



We should expect these to follow simply from the generalized residue theorem just as we 
found for 7 particles. Following the same steps as we did for the two-complex variable case, 
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the generalized residue theorem in the present case is of the form 

J>'i,j 2 ,O = (86) 

i 

Using the explicit form of the box coefficients as residues, the IR equation becomes 

i{l,2,3} + {l,2,4} + i{l,2,5} + {l,2,6} + i{l,2,7} 

+{1, 2, 8} - ^{2, 5, 6} - ^{4, 1, 2} - ±{4, 5, 6} + ±{5, 6, 1} 

+{5, 6, 2} + 1{5, 6, 3} + {5, 6, 4} + ^{5, 6, 7} + {5, 6, 8} 

-i{6, 1,2}- ~{8, 1, 2} - ~{8, 5, 6} = (87) 

Combining terms keeping in mind the antisymmetry of the residues, this follows from com- 
bining the residue theorems as 

^{1,2, i} + ]T{5,6,*} = (88) 

i i 

3.4.2 Two-Loop Evidence 

Not all of the (®) = 56 residues we have are determined by the 1-loop leading singularities. A 
quick look at the map between residues and boxes reveals that all residues which appear at 
1-loop are either of the form {i, j} or {i, i + 3, i+6}, that is, the indices appearing in the 
residues have at least one pair separated by an odd integer. Objects of the form {i, i + 2l, i + 
2m} with all even differences do not appear. Now, this is not quite an invariant statement; 
after all, we can use the residue theorem identities to trade some residues for others. However 
it is easy to see that, beginning with the form of the 1-loop leading singularities we have 
identified, any re-writing would involve sums of pairs of the {i, i + 21, i + 2m} residues. 
Consider an arbitrary residue theorem ^2 k {ji, J2, k} = 0; without loss of generality we can 
put ji = 1. Unless j'2 is also odd, this residue theorem will not contain any of our "missing" 
residues and is irrelevant. If j 2 is odd, then the terms with k odd correspond to our residues, 
and there are even number of these (in fact only two of them since the terms with k = ji, J2 
vanish) . Thus this simple parity argument shows us that we can never express a single residue 
of the form {i, i + 2l,i + 2m} in terms of the others; these objects are not determined by the 
1-loop leading singularities. 

It is therefore natural to conjecture that these residues are associated with a genuine two- 
loop leading singularity. Before making this connection, we can give a natural interpretation 
to the rational functions that appear in these missing residues. As we will discuss in more 
detail in [47], there is a canonical way of beginning with an object and adding a particle 
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to it by applying what we call an "inverse soft factor" [H]. Consider an n particle object 
Ojf ", where we have marked two consecutive colors a, b as being special. We can define an 
(n + 1) particle object O^+i that can be thought of as inserting the particle c between a 
and b. Clearly this operation will have to involve deforming the momenta of a, b in order for 
the new object to conserve momentum with the addition of c. In [17] we will describe this 
operation in a fully supersymmetric way, but for now we will only talk about the addition of 
gluons. We add a positive helicity gluon c + by defining the object 



On+ii-^a, A a ; A c , A c ; Xb, As, • • • ) — ^ ^ 0° b (A a , A' a ; A fc , A' b , ■ ■ ■ ) (89) 

(ac) (co) 



where 



?/ _ (Pa+Pc)\b) , _ (p b + Pc)\b) . . 

{ab) \ba) 

It is easy to see that the deformation preserves momentum conservation with the addition of 
c, and has the correct little group properties as well. In the soft limit p c — > 0, X' ab — > A a 5, and 
O^+i reduces to 0^ b when stripped of the standard soft factor associated with positive 

helicity gluons. This is why we refer to this way of adding a particle as "applying an inverse 
soft factor". A negative helicity gluon is added in the obvious analogous way reversing the 
roles of A and A. 

If we apply the inverse soft factor to the leading singularity associated with a 7 particle 3 
mass box, the resulting object is naturally interpreted as a 2-loop leading singularity, allowing 
us to interpret the missing residues as certain 2-loop leading singularities! The map below 
is associated with adding the particle 8 between 1 and 7 in a 7 particle 3-mass box, which is 
associated with the {3, 5, 7} residue. 




It is obviously beyond the scope of this paper to find all the 2-loop singularities of the 8 
particle amplitude and identify all of them with our residues; as for the lower-point cases we 
have examined, our conjecture makes a clear prediction that all the leading singularities at 
2 loops and beyond are one of the 56 residues we have identified. 
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4 General NMHV Amplitudes 



Using the explicit calculation of the 1-loop box coefficients for all 1-loop NMHV amplitudes 
given in [14J, we have associated all the 1-loop leading singularities with residues. A given 
residue for the n particle NMHV amplitude is determined by giving a list {zi, • • • ,i n -5} 
specifying the minors that vanish. To save space at large n, we can equally well specify the 
minors that are left out and denote the residues by {ji, J2, J3, J4, js}- 

The simplest association of residues with box coefficients is for the 3 mass boxes, which 
are given by a single residue: 




S S-l. 



The "2 mass hard" boxes, with two massive legs and two adjacent massless legs, are also 
very simply identified as the sum of two terms: 




3 2. 



We will refrain from giving the full mapping from residues to the "2 mass easy" and "1 
mass" boxes here, since merely listing the residues is not particularly illuminating. In fact 
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we have found a beautiful structure in this map, which is easiest to study using "Hodges 
diagrams" and other concepts we will more properly introduce in |47]. Elsewhere, we will 
also study the structure of the IR equations more systematically, and show that the same 
pattern we saw for n = 7, 8 persists to all n: the 1-loop IR equations follow from residue 
theorems. 

4.1 Tree Amplitudes 

Here we give a compact discussion of all tree NMHV amplitudes. Since we will be dealing 
with sums over many residues, it is convenient to introduce some algebraic notation for 
doing this. We can denote a given residue {ii,Z2, •'* j^n-s} by an antisymmetric product 
{*i}{*2} • ■ • {in-§}- In this way we can talk about an expression like ({l} + {3})({4} + {7}) = 
{1, 4} + {1, 7} + {3, 4} + {3, 7}. In order to write the NMHV amplitudes in a compact way, 
we introduce a little more notation. Let us define 

{{ii}{i 2 } if k < 12 ) 
\ (91) 
otherwise J 

Let us also define the formal sums "£" (for "even" ) and "0" (for "odd" ) 

£= £{*},£>= (92) 

k even k odd 

The BCFW form of the NMHV amplitude is then given by the following sum of residues: 

in — 5) factors 

Just to be explicit, we can write this out for n = 8 as the 10 terms 

{2}{3}{4} + {2}{3}{6} + {2}{3}{8} + {2}{5}{6} + {2}{5}{8} + {2}{7}{8} 

+{4}{5}{6} + {4}{5}{8} + {4}{7}{8} + {6}{7}{8} (94) 

The P(BCFW) form of the amplitude is of exactly the same form, but exchanging £ «-> —O. 
Since the cyclic shift also exchanges £ and O, the expected identity M B cfw — ^p(bcfw) will 
also guarantee the cyclic invariance of the amplitude. We thus turn to proving this identity. 

4.2 Cyclic Invariance and Homology Classes 

Let's start by writing the general form of the residue theorems, which in a completely straight- 
forward generalization of what we have seen already is of the form 

D»iK«3}-{UW = ( 95 ) 
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We will show that these identities imply the desired statement 



£*£>*•-• 



(l) n - 5 (D*£ 



(96) 



Note that the residues involved in the left and right hand sides are completely distinct, and 
the unphysical poles on the left and right hand sides are different. Thus aside from cyclic 
invariance, this identity also enforces the absence of unphysical poles as well. 

Before giving the proof in the general case, let us go back to the case of n = 7 for 
simplicity, and ask the general question: under what conditions is the object 



cyclically invariant? Most naively, cyclic invariance requires /j+i^+i = however, this is 
too strong a requirement. The reason is that because of the residue theorems, we can add to 
F the quantity zero in the form 



Thus / is ambiguous; because of the identities two sets of / give the same object F: 



This highlights an important point. As we have seen repeatedly, the amplitude is not associ- 
ated with a unique contour of integration picking out a unique set of residues, but rather by 
an equivalence class of residues defined up to combinations that vanish due to residue theo- 
rems. We can say this more intuitively and geometrically by saying there is an equivalence 
class of contours that can smoothly be deformed into each other; in this sense the amplitude 
is associated with some homology class. Now the homology class for the tree amplitude 
should be cyclic, even though any individual representative may not be manifestly cyclic. 
It would be clearly be nice to be able to understand the cyclic symmetry properties of the 
homology class directly. 

In order to do this, let's interpret the redundancy of equation ( 199]) as a sort of "gauge 
invariance"; for the object F to be cyclically invariant it is not necessary for / to be cyclically 
invariant, only that it is cyclically invariant up to a gauge transformation. As usual, it is 
preferable to deal with "gauge-invariant" quantities. In this case, if we define the difference 
operators Ai >2 via 




(97) 




(98) 



fj,k — > fj,k + — Oik 



(99) 



^■lfj,k — fj,k ~ fj-l,ki A 2 fj t k — fj,k — fj,k-l 



(100) 
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then it is obvious that the "gauge invariant" object made out of / is 

(AxAa)/ (101) 

This "curvature" is an invariant characterization of the homology class defining F. And in 
particular, the quantity F is cyclically invariant only if the curvature AiA 2 / is cyclically 
invariant; the converse is also easy to prove so this is an if and only if statement. 
Now it is trivial to see that 

(AxA 2 ) {a} * {b} = [{a} - {a - 1}] * [{b} - {b - 1}] (102) 

So we can look at the curvature associated with the BCFW form of the amplitude: 

(AxA 2 ) [({2} + {4} + {6}) * ({3} + {5} + {7})] = 

[{2} - {1} + {4} - {3} + {6} - {5}] * [{3} - {2} + {5} - {4} + {7} - {6}] 
= - [{1} - {2} + {3} - {4} + {5} - {6} + {7}] * [{1} - {2} + {3} - {4} + {5} - {6} + {7}] 

where in the last line we added (7) to the first bracket, since (7)*(a) = 0, and also (1) to the 
second bracket, since (a)*(l) = 0. The last line is then manifestly cyclic invariant, since we 
have a single factor "squared", and this factor goes into minus itself under a cyclic shift. 
This pattern continues trivially for all n; the "gauge invariant curvature" is 

(A 1 ...A n _ 5 )f ju ... Jn _ 5 (103) 

The "curvatures" associated with the BCFW and P(BCFW) forms of the amplitude can then 
easily be seen to be identical 

(Ax-.-A^s) [5*0*...] = (-l) n - 5 (A!...A n _ 5 ) [0*£*---] = -[(£-0)*(£-0)*---] 
which establishes that Mbcfw = -^p(bcfw)- 

5 Residues in Multidimensional Complex Analysis and 
the Global Residue Theorem 

In moving beyond NMHV amplitudes, with k > 3, we encounter a number of new features. 
Each minor becomes a higher than linear degree polynomial in several r variables. We 
are interested in making sense of residues and residue theorems in this general situation. 
This is a basic subject in multi-variable complex analysis and algebraic geometry. The 
identification of residues is entirely analogous to the simple treatment we gave in discussing 
NMHV amplitudes. The procedure for deriving residue theorems by repeated application of 
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the ordinary Cauchy theorem one variable at a time does not extend to the general case; 
however there is still a residue theorem of exactly the same type we have seen, known as the 
global residue theorem. The subject of residues in multidimensional complex analysis has 
not made very many appearances in the physics literature (the only example we are aware 
of is [49J). For this reason in this section we provide a short, self-contained review of the 
subject in the generality needed for the purposes of this paper. Most of the discussion will 
follow EH and EH. 



5.1 Local Residues 

Consider a holomorphic mapping / = (ft, ■ ■ ■ , f n ) '■ C n — > C n . For us each f i is in fact a 
polynomial of a fixed degree. Let us assume that / has an isolated zero at a 6 C n . In other 
words, if we choose a small enough neighborhood U of a, then / _1 (0) D U = {a}. The local 
residue of a meromorphic form lu — h{z)dzj f\ . . . f n with dz = dz\ A ... A dz n at the point a 
is defined by the integral 



over the contour 



T a = {zeC n : \fi(z)\ = ej, j = 1, . . . , n}, (105) 



h(a) 

resMa = with J f (z) 

Jf(a) 



where €j G 1R + are small positive real numbers. Note that this is the straightforward gener- 
alization of the natural construction we used for 7 particles^. 
The value of this integral is given by 

9(fl, ■ ■ ■ , fn) ClDfil 

d(zi,...,z n ) 

The most important piece of information we did not specify carefully in our 7-particle 
example was the definition of the orientation of the contour. Now that we have introduced 
the form notation we can define the orientation unambiguously to be such that d{axgfi) A 
. . . A d(axgf n ) > 0. More precisely, the cycle under consideration can be parameterized as 
fi(z) = eje 1 ^, with 8i = arg/j. Therefore, the orientation is determined by the order of the 
#j's in the measure. In a sense one can think of (11041) as the local residue of h(z) with respect 
to the mapping /. It is in this sense that the local residue is skew-symmetric with respect 
to permutations p of the components of /: 



res 



(/ lr ../„) = (sgn p) res (/p(1)v .. Jp(n)) . (107) 



3 Note that n here denotes the dimension of the space of free integration variables and not the number of 
particles! 
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As mentioned in the seven-particle case, the (S 1 )™ ^ T n contour does not enclose the 
point a simply because it has real dimension n and one needs a 2n — 1 dimensional subspace 
to do that (only for n = 1 the two dimensions agree). It is interesting to note that there is 
an equivalent definition of the residue in terms of an integration over a S 2 "^ 1 sphere that 
encloses the point a. The integral is of a 2n — 1 real form defined in terms of / and /. We 
will not make use of this and we refer the reader to [50J or [51J for details. The reason we 
mentioned this fact is that it enters in the proof of the following results which we simply 
state. The proofs follow intuitively by using the S 2n ~ l spheres and Stokes' theorem. 

5.1.1 Total Sum of Residues 

As we discovered in the seven-particle discussion there is a natural generalization of the 
residue theorem to the multidimensional case. Let us state the theorem in the generality we 
need. 

Let uj = hdzj fx .../„ be defined by polynomials h and /j. Let Fj = {z e C n : fi(z) = 0} 
be the hypersurface (i.e. n — 1 dimensional subspace) associated with and Z = F± n F 2 fl 
. . . fl F n be the set of zeroes of /. Here we assume that Z is a discrete set of points. Then 
one defines the Global residue of h with respect to the map / as 



R eS/ (h) = ^res(cu) a . (108 



Now, the Global Residue Theorem (GRT) states that if deg(h) < deg(/i) + . . . + deg(/ n ) — n 
then 

Res f (h) = 0. (109) 

Just as we found in section 3.3.4, this theorem allows us to find relations among the 
local residues of a form uj with / : C n — > C m and m > n. Once again we assume that 
Z a = F ai fl . . . fl F an is discrete for any subset a = {cui, . . . a n } C {1, . . . , m}. We also 
assume that Z a fl Zp = for a ^ (3. For each Z a we can define 

Res Q uj = res(u;) a . (HO) 

aeZ Q 

Note that the sum of the residues Res a (u;) is skew-symmetric with respect to «i, . . . , a n . 
Now, for each set 7 = {71, . . . , 7„_i} C {1, . . . , m} one has 

Res htj) u = (111) 

ie{l,...,[7],...,m} 

where [7] in the range of the sum means excluding the set 7. This follows easily from the 
GRT. 
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A consequence of this is that for any partition of the set {1, . . . , m} into n disjoint subsets 
Ji, . . . , J n we also have the relation 

E Res (n,-^ = (112) 
jieJi,...,j„eJ„ 

We can summarize the implications of the global residue theorem for our particular ap- 
plication as follows. Let D = (k — 2) (n — k — 2) be the number of free variables. To define a 
residue, we solve the equations putting D minors (ji • • • ji+k — 1) = ■ ■ ■ = (jo • • • jD+k — 1) = 
0. Since these are in general polynomial equations of degree higher than one, there will be 
many solutions, which we index by a variable A. The corresponding residue is denoted as 
{ji, ' ' ' j 3d} a- The global residue theorem is then the statement that for any partition of 
the set {1, • • • , n} into D disjoint subsets Ji, • • • , Jd, 

E E <•>'• ./"!.* ° ( 113 ) 

A jieJi,...,j n e.JD 

as long as there are no poles at infinity. Due to the antisymmetry property of the residues, 
these all follow from taking linear combinations of the set of equations 

EE^---'^-i>o = o (H4) 

A i 

Note that in this discussion we have tacitly assumed that n > D in order to even be able 
to define the residues. It is however clear that for large enough n ~ 2k that D ~ k 2 ^> n. 
We will discuss how to think about residues in this case in section 7. 

6 N 2 MHV Amplitudes 

Up until now we have only studied NMHV amplitudes. There are two features of NMHV 
amplitudes which make them special. The first is that for k = 3 each factor in the de- 
nominator is linear. The second is that at one-loop it is easy to see from the quadruple 
cut formula for box coefficients that all four-mass boxes vanish trivially. In this section we 
will study N 2 MHV amplitudes, which are the simplest case where the minors are not linear 
but quadratic polynomials in 2(n — 6) variables. Most of our discussion will focus on the 
8 particle amplitude, which is also the simplest case in which four-mass box coefficients are 
non-vanishing. At the end of the section we will also take a peek at the 9 particle amplitude. 

This section is divided into five parts. In the first we concentrate on showing how a 
four-mass box coefficient appears as a residue in our formula for the 8 particle N 2 MHV 
amplitude. Four-mass boxes are special because they are completely IR finite; furthermore 
unlike all other boxes we have seen so far, where each of the two leading singularities have 
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been rational functions of the kinematical variables, for the 4 mass boxes each leading singu- 
larity has a square-root dependence on the kinematical variables. Even more than our explicit 
identification of all NMHV box coefficients, this gives smoking-gun evidence that our conjec- 
ture is correctly computing 1-loop leading singularities. Since our main goal here is simply 
to verify the presence of 4 mass boxes, in this part we study the simplest "split-helicity" 
configuration. 

In the second part, we concentrate on the alternating helicity configuration and map out 
all residues coming from setting four of the eight denominators to zero. In other words, we 
solve (T\ = 70 equations; we find that each set of equations has two solutions, giving a total 
of 140 residues. Many of the objects familiar from the tree and 1-loop alternating helicity 
amplitudes can be identified with these residues; in particular we identify the objects T, U, V 
introduced in [2] to express the BCFW form tree-level 8-point amplitude, as well as P(T),P({7) 
and P(V) in the P(BCFW) form [17J. Using the same "curvature" notions we used in our 
analysis of general NMHV tree amplitudes, we can easily derive the Mbcfw = A^p(bcfw) 
equality, which is a remarkable 40 term identity. 

Moving on, however, we find that some of the one-loop leading singularities are no-where 
to be found in this list of 140 residue. These include the two other objects W, X that naturally 
occur [2] in the BCFW form of the tree-amplitude. Where could these objects be hiding? 
In the third part of this section we will see the answer is that they are associated with a 
new kind of residue. Mathematicians refer to these as "composite residues", and indeed 
their appearance will be ubiquitous for all k > 4 and n > 10; in this section we will give 
a very simple intuitive explanation of what these residues are and how to find them, and 
defer a general discussion of their properties to the next section where the residues of maps 
/ : C n — > C p with p < n are considered. Armed with the more general notion of composite 
residues, we proceed to identify all the 1-loop box coefficients for the 8 particle N 2 MHV 
amplitude. 

However, just as we found for the 8 particle NMHV amplitude, in the fourth part of this 
section we show that even some of the "usual" residues are not associated with any 1-loop 
leading singularities. It is again natural to think that these are naturally computing 2-loop 
leading singularities. There is one more piece of evidence in favor of this interpretation. Up 
to now we have seen that the residue theorems are associated with 1-loop IR equations. Now, 
there are clearly global residue theorems involving the residues associated with 4 mass boxes. 
However since the 4 mass scalar boxes are completely IR finite, they should not make any 
appearance in 1-loop IR equations. Interestingly, however, we find that precisely the residue 
theorems involving the 4 mass boxes also involve the "mystery" residues with no 1-loop 
interpretation. The natural interpretation is that these are in fact two loop IR equations, 
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which relate the IR divergent part of the two-loop amplitude to the 1-loop amplitude! An 
amusing application of this residue theorem is an expression for the 4-mass box coefficient 
which is explicitly a rational function of the kinematical invariants. 

We end by taking a brief look at the 9 particle amplitude, and by showing that the set of 
residues include information that is unambiguously absent at 1-loop, further reinforcing our 
claim that all leading singularities at all loop order are being correctly computed. 

6.1 Four-Mass Box Coefficient 

Let us start with a brief comment about the computation of the coefficient using the quadruple 
cut approach. Consider the four-mass box with external legs given by K\ = p$ + pi, K 2 = 
V2 + P3, K 3 = pi + p 5 and K4 = p e + p 7 . Cutting all four propagators of the box give rise to 
the following equations for the loop momentum £, 

£ 2 = 0, (t - K x f = 0, (I - K x - K 2 f = 0, (£ + K 4 ) 2 = 0. (115) 

One way to solve these equation is to parameterize £ a a = (A3 + «A4) a (/9A3 + 7A4)q. Plugging 
this into the equations we find two solutions obtained from solving a quadratic equation. The 
discriminant of the quadratic equation is given by 

K 2 K 2 K 2 K 2 
A = 1 - 2( Pl + p 2 ) + ( Pl - p 2 ) 2 with Pl = — i-f- and p 2 = (116) 

-^12-^23 iV 12- fV 23 



This means that a, (3 and 7 are of the form a ± by A for some a and b rational functions of 
the kinematical invariants. 

For simplicity consider the one-loop amplitude M(l~, 2~, 3 _ , 4~, 5 + , 6 + , 7 + , 8 + ). In this 
case the box we are studying is the only non-zero one. A non-trivial test of our proposal 
is that out of all residues of the integral with denominator (2345) (3456) . . . (8123) only one 
combination of four factors should give rise to residues located at points with a square root 
dependence on the kinematical invariants. Indeed we have checked that (2345) = (4567) = 
(6781) = (8123) = is the only combination that gives a quadratic equation and therefore 
two residues located at points that depend on a square root. The argument of the square 
root turns out to be precisely equal to A! 

The coefficient we are interested in is computed in the quadruple cut approach as the 
product of four four-particle tree-amplitudes (normalized with the Jacobian of the 4-mass 
box). As a sample let us write two of the four amplitudes 

(23> 3 (^-P2 3 4) 3 



(3 t - p 23 ) (i - P23 t) (i 2) ' (45) (5 t - p 2U5 )(£ - p 23 45 i ~ P23) 



(117) 
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The product of all four amplitudes can be simplified to give B(£) as 

(23) 3 [67] 3 [4|l-p 23 |4) 3 (ll) 3 

<45>(81>(£ 2)<3|£-p 2 3|4][4|p 2 3|£><€|p 18 |7]<5|^-p 2 345|6]<8|(£-p 18 )(^-p 23 45)(^ -P 23 )|4] 

The coefficient of the box is then B(£^) + B{£S~^) where £ ± are the two solutions to the cut 
equations. 

We have numerically verified that the residues at the two points in C 4 which are solutions 
of (2345) = (4567) = (6781) = (8123) = match B(£^) and B{£^) respectively! Extending 
our now familiar notation for residues, we denote these two solutions by {2468}i and {2468}2, 
and so 

£(£(+)) = {2468}a, B(£^) = {2468} 2 (118) 

A final comment is in order here. The global residue theorem implies that B{£ <y+ ' ) )+B{£ <y ^) 
is equal to the sum of other residues. As an illustration let us choose fi = (2345), / 2 = (4567), 
f 3 = (6781) and / 4 = (8123) (3456) (5678) (7812). Recall that (1234) = 1 in the natural gauge 
choice for this helicity configuration. The residue theorem then gives 

B{£^) + B{£^) = -J2 ({2463} A + {2465} A + {2467} A ) (119) 

A 

As we have mentioned, all the residues other than the ones determining B{£^) are rational 
functions of the kinematical invariants; they are determined by solving linear equations and 
do not contain square roots. We therefore obtain a remarkable form of the 4 mass box 
coefficients directly as a rational function of the kinematical invariants. 

6.2 Tree Amplitude 

We now turn our attention to discussing the alternating helicity amplitude with assignment 
l + 2~3 + 4~5 + 6~7 + 8~, beginning with a discussion of the tree amplitude. In [2], the 20 terms 
of BCFW form of this amplitude were determined in terms of five basic objects T, U, V, W, X, 
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together with their cyclic images. Explicitly, these 5 objects are given by 

rp = [13] 4 (46) 4 (68) 4 

[12] [23] <45> <56> (67) <78) <6 14+5 13] <6 1 7+8 1 1] <6| (7+8) (1+2+3) |4> <8 1 (1+2+3) (4+5) |6> 
v _ [13] 4 [57] 4 (48) 4 



V 



[12][23][56][67](4|2+3|l](4|5+6|7](8|l+2|3](8|6+7|5]tf ] 4 3] 

[13] 4 (46) 4 (8|l+2+3|7] 4 

[12][23](45)(56)[l|(2+3)(4+5+6)|7](4|5+6|7](8|l+2|3](8|(l+2+3)(4+5)|6)tf ] 4 3] 4 4] 

, 1 , 0|r,l4 /or,\4 



_ [35] 4 <6|8+1+2|3] 4 <82> 



[34] [45]<12>(67>(81)[5|^ 2] i^[ 3] |3]{2|i^| 3] (6+7)(4+5)|3]<6|4+5|3]<7|i^| 3] |3]<8|l+2|^ 



x _ [35] 4 (82) 



[34] [45] [56](78)(81)(12)(2|3+4+5|6](7|8+l+2|3]4 4] 

(120) 



2 



Here we use the notation of [2] where = Pi~\ — ■ +Pi +m -i and r™ 11 = (p$ + - • ■ +Pi +m -i 
In [2] it was also noted that, using a peculiar- looking 6 term identity, the tree amplitude 
could be written purely in terms of T, U, V, together with their cyclic images. As we will 
explain in [17], this identity is in fact nothing other than our 6 term identity for the 6 particle 
amplitude, dressed into a 6 term identity for the 8 particles by the application of two "inverse 
soft factors" . At any rate, the nicest form of the alternating helicity amplitude is then given 
by 

M^ w + - + - = (l + g + --- + g 7 )[T + V] + (l + g + --- + g 3 )U (121) 

Here g is the operation that sends i — > i + 1 and conjugates ()•*->•[]. We have only 4 images 
of U since it has a flip symmetry. The P(BCFW) form of the amplitude is of exactly the 
same form, with T, U, V ->P(T),P(f7),P(V), given by 

p{g , T) ([75](2|7+8|6] + [76](2|1+3|5]) 4 



P(g 4 u) 
P(g 4 v) 



[45] [56] [67] [78] < 12) <23> <3 1 4+5 1 6] < 1 1 8+7| 6] [4| (3+2+ 1 ) (8+7) 1 6] [6 1 (5+4) (3+2+1 ) | i 

(2|1+3|5+7|6) 4 



(12)(23)(56)(67)(l|2+3|4](7|6+5|4](3|2+l|8](5|6+7|8]tf ] 4 3] 

(2|(l+3)(8+l+2+3)(4+6)|5] 4 

[45][56](12)(23)(7|(6+5+4)(3+2)|l)(7|6+5|4](3|2+l|8][6|(5+4)(3+2+l)|8]4 3l 4 3l 4 4] 

Clearly, the expected identity 

JW BCFW — JW P(BCFW) 

is an amazingly non-trivial 40 term identity! 
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Let us try and identify these with residues. As we mentioned before, the residues are 
associated with putting 4 of the minors to 0; this amounts to solving four coupled quadratic 
equations. Somewhat surprisingly, it turns out there are always two solutions to these equa- 
tions. We will understand this fact more deeply when we discuss the relationship of our 
conjecture to algebraic geometry and the Schubert calculus in section [HJ This means that 
we should label the residues as {zi, z 2 , 13, 24)1,2, where the subscript denotes the two solutions 
associated with putting the four minors beginning with z 1( z 2 , is, to zero. 

Now, we find a very pleasant surprise in identifying T, U, V and P(T),P(i7),P(V) with 
residues. A term and its Parity conjugate are the two different solutions associated with a 
single residue! Specifically, we find 

T = {3,4,5,6}! , P(T) = -{3,4,5,6} 2 
U = {4,5,8,1} 1 , P(f/) = -{4,5,8,l} 2 

^ = {4,5,8,3}! , P(y) = -{4,5,8,3} 2 (123) 

Under the conjugation operation g, there is a minus sign in the identification, so that e.g. 
gT = —{2, 3, 4, 5}. The 40 term identity then becomes 

£ K 1 - 9 + ■ ■ ■ - 9 7 ) ({3, 4, 5, Q} A + {4, 5, 8, 3} A ) + (1 -<? + -••- £ 3 ){4, 5, 8, 1} J = 

A=l,2 

This looks exactly like a residue theorem! In fact, the global residue theorems for our case 
are linear combinations of the statements 

J2 ihJ2,j3,i} A = (124) 

i;A=l,2 

Thus in order to check whether a sum of the form 

^ = EEWM^,0i (125) 

A i,j,k,l 

vanishes, it suffices to see whether the "curvature" vanishes 

(A 1 A 2 A 3 A 4 )/ iJ - = (126) 

A completely straightforward computation shows that the "curvature" of the summand in 
the 40 term identity indeed vanishes, establishing the validity of M B cfw = ^p(bcfw) m this 
case. 



48 



6.3 A New Class of Residues 



We have identified T, U, V and their parity conjugates, which are everything needed for the 
tree amplitude. However, W, X do appear directly in the BCFW formula, and are some of 
the 1-loop leading singularities. Indeed W, its flipped cousin rW (where r : % — * 9 — %) and 
their parity conjugates are given by 



W = {7, 8, 4, 6}x, P(W) = {7, 8, 4, 6} 2 ; rW = {3, 4, 6, 2} u P(rW) = {3, 4, 6, 2} 2 (127) 



while the parity conjugates are given by the second roots. However, examining all 140 
residues, we could not find any objects that matched X\ Where could X possibly be hiding? 

The answer is quite remarkable, and illustrates the rich structure of our conjecture. The 
definition we have given for residues is valid for completely generic rational functions of many 
variables. However we have a very special set of functions, and an interesting phenomenon 
can occur that allows us to define residues in a more interesting way. We can illustrate the 
main idea by considering a function of 3 complex variables x, y, z of the form 



There are only two polynomial factors in the denominator while the definition of residue we 
have been using would require three polynomials to determine the three variables. However, 
something special happens on the locus where the first polynomial x vanishes: the second 
polynomial (x + yz) — > yz factorizes, so it seems sensible to define the residue of this function 
at x = y = z = to be 1. 

We find exactly the same structure in our problem! We find that in some cases, in putting 
only two of the minors (i, i + 1, i + 2, i + 3) = and (J, j + 1, j + 2, j + 3) = 0, a third minor 
(k, k + 1, k + 2, k + 3) factorizes into a product of two linear factors in the four r's. Therefore 
we can define a residue using only these three minors. These residues are called "composite 
residues" by mathematicians and will be discussed more systematically in the next section. 
We find that there are still two composite residues. We use a notation {i,j, k 2 }i^ 2 to denote 
these residues. 

In order to see the composite residues explicitly in our example, it is easiest to use the 
gauge fixing where 



1 



(128) 



x(x + yz) 



C = 



I 1 C15 Cie C17 Cis \ 

1 c 25 c 26 c 27 c 28 

1 c 35 c 36 c 37 c 38 

\ 1 c 45 c 46 c 47 c 48 / 



(129) 



49 



Consider the first four determinants of the minors that enter in our formula, 



(1234) 
(2345) 

(3456) 



(4567) 



1, 

c 15; 

Cl5 Ci6 
C25 C 26 



Cl5 Ci 6 C17 
C25 C 26 C 27 
C35 C 36 C37 



(130) 



We will now show that, it is possible to define a residue using only three minors (2345), 
(3456), (4567). Recall that the en are linear in the ta- Now, note that on the locus of 
(2345) = 0, (3456) factorizes into the product of two terms: 



(2345) = 



C15 = 



(3456) = -ci 6 c 25 



'131} 



So just using these two minors allows us to solve for 3 r's, by setting C15 = c\q = C25 = 0. Then 
(4567) is a quadratic polynomial in the remaining r variable, which we can set to zero finding 
two solutions. The residue is defined in the usual way, with the relevant Jacobian being the 
one determined by the four objects we are setting to zero <9(ci 5 , Ci 6 , c 2 5, (4567))/<9(ti, • • ■ , r 4 ). 
The two composite residues in this example are denoted by {23 2 4} lj2 . 

Having identified all the sorts of residues we should be looking for, we can now identify 
all the 1-loop leading singularities! For instance the object X is 



X = {6,7 2 ,8} ] 



(132) 



We can now give the map from residues to all 1-loop leading singularities. We have implicitly 
given all the coefficients for the 2 mass hard and 1 mass boxes in giving the BCFW and 
P(BCFW) form of the tree amplitude. The four mass boxes are given as 



7. 
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The 3 mass and 2 mass easy boxes are given as 



7. 



7 



1 



T 



2 t 



cxVW 



7 



1 



9_ 



7. 



1 



^3 

5 "f 



where here the operation r is a reflection, sending i — > 9 



z. 



6.4 Two Loop Leading Singularities 

Just as we saw in the 8 particle NMHV amplitude, even amongst the 140 "usual" residues, 
there are objects we can not identify with any 1-loop leading singularities, of the form 
{1,2,4, 6}a and {1,2, 5, 7}^ together with their cyclic images. It is again natural to as- 
sociate these with two-loop leading singularities. Following the same strategy as for the 
NMHV case, we look for these residues by applying "inverse soft factors" on the 7 particle 
1-loop leading singularities, and then associate these with a certain 2-loop leading singularity. 
Doing this we find the mapping 
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There is one more piece of evidence that these residues are associated with two-loop 
physics. Let us return to the global residue theorem involving the 4 mass box coefficient 
{1, 3, 5, 7} a, which we can write as ^. A {1, k, 5, 7} = 0. Let us group the terms as 

3, 5, 7} A + {1,6, 5, 7} A + {1,8, 5, 7}^) = ^ (-{1, 2, 5, 7} A + {4, 5, 7, 1} A ) (133) 

A A 

The second two terms on the LHS are just cyclic images of W, so the left hand side groups 
together objects that appear at 1-loop, while the RHS has the residues with no known 1- 
loop interpretation. What could the physical interpretation of this residue theorem be? So 
far, all residue theorems have been related to IR equations, but as we have stressed the 
4 mass box is completely IR finite and would not participate in any 1-loop IR equation. 
However, recognizing the terms on the RHS as giving two-loop singularities allows a very 
natural interpretation: the residue theorem should be thought of as a 2-loop IR equation, 
which relates the 2-loop IR divergences to the 1-loop amplitude. IR finite 4 mass boxes can 
then occur as a part of this 1 loop amplitude, just as the 1-loop IR equations relate the 
IR divergent part of the 1-loop amplitude to the IR finite tree amplitude. It would be nice 
to check this picture in detail. The 2-loop IR equations have been understood by Catani 
[53] ; proceeding further would require us to identify all the 2-loop leading singularities and 
reconstruct the amplitude from them, which is beyond the scope of the present work. 

Finally, as we already mentioned, this global residue theorem allows us to write the 4- 
mass box coefficient in a new way, manifestly as a sum of rational terms. We will not give the 
explicit expression here, and simply mention that all the terms appearing in this equation are 
explicitly known: the 1-loop residues are related to W and P(W), and the 2-loop residues are 
obtained by applying an inverse soft factor, defined in section 3.4.2, to the 7-particle 3-mass 
box coefficients. 

6.5 9 Particle N 2 MHV amplitude 

We end with very brief comments about the 9 particle amplitude. We will not present 
a detailed mapping between residues and leading singularities; instead our goal is to give 
another sharp proof that the residues contain information beyond one loop. 

The strategy for doing this is simple. The best way of describing it is to let all the initial 
data encoded in the spinors A and A of each particle to be rational numbers, i.e, to think 
about the n- vectors defining the A and A-planes as vectors in Q". Suppose that all the 
residues are secretly associated only with one-loop leading singularities. Now, we know that 
the one-loop leading singularity is a product of tree amplitudes, which are rational functions 
of the kinematical invariants. However the values of the cut momenta £ can be irrational; 
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indeed as we remarked in our discussion of the 4 mass box coefficients for the 8 particle 
N 2 MHV amplitude, £ ± = A ± B\fK where the discriminant A is reproduced below: 



A = 1 - 2( Pl + p 2 ) + ( Pl - p 2 ) 2 with Pl = and p 2 = (134) 

-^12^23 -^12-^23 

Note that pi and p 2 belong to Q by assumption. 

This says that for the residues associated with 1-loop leading singularities, the number 
field E of solutions for the residues can only be a quadratic extension of Q, moreover, the 
only square roots that appear can only be of the form v^A for some box. Furthermore, with 
9 particles, there are only 9 four-mass boxes and these are cyclically related, so the list of 
possible square roots that can appear is a small one. We have checked all the usual (non- 
composite) residues, and found that only {1,3,5,7,8,9}^ and its cyclic partners contain a 
square root, and that the argument of the square root is precisely the A associated with the 
box. However moving to the composite residues, we find a wide array of "new" square roots, 
which can definitely not be expressed in the number field only containing rationals and the 
square roots of the one loop A's. This proves conclusively that the residues contains a great 
deal of information beyond one loop. 



7 Composite Residues and Large k Amplitudes 

As we already mentioned at the end of section 5, the number of free variables of integration 
in our formula, (k — 2) (n — k — 2), is in general larger than the number of determinant 
factors in the denominator which is n. This means that for large n and k, the "ordinary" 
definition of residues we have been using, which requires at least as many polynomial factors 
as complex variables, can not be used. Instead, the new sort of "composite" residues we 
already encountered in the previous section will be needed. That is, we will find that on the 
zero locus of a given minor, the subsequent minor factors into pieces allowing us to solve for 
more of the r variables. In the first part of this section we show that doing this allows us to 
solve for all the r variables and determine residues. We leave a more complete analysis, as 
well as a physical interpretation of these residues, to future work. In the second part of this 
section, we briefly discuss the way in which these residues are discussed in the mathematical 
literature. 
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7.1 The Large n = 2k Limit 



Let us concentrate on the most extreme case, i.e., consider n = 2k. In this case we have 
(k — 2) 2 integration variables. Consider once again the gauge fixing where 

• Ci ; fe + i • • • c l Ci 9k- \ 

■ C 2 ,fc+l C 2 ,fc+2 • 



c = 



/ 1 
1 







1 Cfc,fc+1 Ck,k+2 



Cl,2fe-1 
C2,2fc-1 



Cl,2fc 
C2,2fc 



Cfc,2fc-1 Cfc, 2 fc / 



(135) 



Following the same analysis as in the eight particle case, let us list the first four of the 
determinant factors, 

(12 ... A;) = 1, 
(23... k + 1) = ci ifc+ i, 



(34 



A; + 2) = 




Cl,fc+2 






C2,fe+1 


C 2 ,A:+2 




A; + 3) = 


c l,fc+l 


C\,k+2 


c l,fc+3 


C2,fe+1 


C2,fe+2 


C2,fc+3 




C3,fe+1 


C3,fe+2 


C3,fc+3 



(136) 



Let us denote the free variables by r = (ti, . . . , £ C^ -2 -* 2 . Each factor is at most a 

polynomial of degree (k — 2) 2 — 2. Now, let's see what happens as we begin to set minor 
factors to 0. As we saw already, 



(234 ... (fc + l)) =0 (345 . . . (k + 2)) = -c lik+2 c 2 , k+1 



(137) 



Note that ci jk +2C2,k+i is the product of the elements on the "secondary" or "skew" diagonal 
of the minor (345 ... (k + 2)). Now, putting all three of these factors to zero, we find 



(456 ... (k + 3)) = c ltk+3 c 2 ,k+2C3,k+i 



(138) 



Now it is easy to see the pattern. Once we use all factors in a given step, the determinant of 
the next minor reduces to the product of the entries in the secondary diagonal of the minor. 

Note that the matrices increase in size until one of them becomes a k x k matrix and then 
start decreasing. At this point one has to be careful to remember that even though we seem 
to get k linear factors coming from the k C/j's in the secondary diagonal, we know from the 
general analysis in section 4.3 that the order cannot be larger than k — 2. This means that 
2 of the factors become trivial, i.e., r-independent. 

Now we can count the number of factors we can set to zero to define residues. Using the 
first k — 2 determinants (excluding (123 ... A;)) we find 1 + 2 + . . . + (k - 2) = (k - l)(k - 2)/2 
factors. Using the last k — 2 determinants we also find the same number. Therefore we 
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have (k — l)(k — 2) factors which is larger that then number of free variables. This gives 
many different residues. Once again, the precise counting of these residues as well as their 
precise physical interpretation falls outside the scope of this paper and it is left for future 
investigations. 

7.2 Composite Residues and Residual Forms 

Mathematicians naturally encounter composite residues in studying mappings / : C n — > C p 
consisting of p irreducible polynomials in n variables, with n > p. It is clear the previous 
discussion of residues does apply straightforwardly in this case. In pjTJ the concept of residual 
currents and in [52J residue forms are introduced as a way to study these cases. We find the 
residue forms to be naturally adapted for our problem. Let us give the definition of a residue 
form. Let S = {z G C n : s(z) = 0} and u)(z) = (h(z) / s(z))dz with h(z) and s(z) holomorphic 
and dz = dz\ A ... A dz n , then the residue form is an (n — l)-fornQ given by 



tf-i 




dz m (139) 

s 



at the points where s' z , ^ 0. Let us explain the different objects entering this formula. The 
formula is defined for any j such that the condition is satisfied, dzy] is a n — 1 form obtained 
by wedging all dz^s except for dzj. The function h{z)/s' z . is evaluated on S(z) = and 
therefore it depends on only n — 1 variables. Finally, s' z . is the partial derivative with respect 

to Zj. 

Now one can define the (n — m)-form 

res m [o;] = res m o ••• oresi[c<j] (140) 

called the composite residue form. 

In simple cases, like when one has a map / : C n — > C n with isolated zeroes, the composite 
residue form res n [a;] coincides with res(a;) in (11041) and different orders of computing the 
composite residue form give rise to the same answer up a sign, but this is not true in general. 
Using this notation, we can return to our simple example of the function Zl ( Zl +^ 2Z3 ) > putting 
s(z) = z\{z\ + Z2Z3) and h(z) = 1. Then, res 2l (co>) = l/z 2 z 3 . The following composite residue 
form gives 

res 3 [a;] = res Z3 o res 22 o res Zl [u] = 1. (141) 

Mathematicians have also identified residue theorems involving composite residue forms. Our 
real interest in making a connection with the mathematical literature here is to exploit these 



*More precisely, we have a (n — l,0)-form. 
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results, since as we have seen residue theorems have some deep physical content in our picture. 
We will take up this subject at greater length in future work. 

We find it intriguing that the notion of composite residue has made two appearances in 
our story so far, first in defining certain kinds of leading singularities as residues in the loop 
integrals, and now, in identifying certain leading singularities as residues in the integral over 
the r 7 variables. Needless to say it would be fascinating if this were more than a co-incidental 
fact! 

8 Amplitudes and the Schubert Calculus 

We have repeatedly mentioned that the pieces entering in our formula all have natural ge- 
ometric interpretations in terms of Grassmannians. However, we have not made use of this 
fact in any of our computations. It is clear that a deeper connection with the algebraic 
geometry of Grassmannians could shed new light into our formula. In this section we give 
a simple example of the use of the homology and its ring structure in order to compute the 
number of solutions to a given set of equations coming from distinct minor factors. 

The first step is to define the basis of homology for a Grassmannian G(k,n). The ring 
structure is obtained by using as multiplication the intersection of cycles. The cycles in the 
basis are called Schubert cycles and the ring structure gives rise to the Schubert calculus. 

In this discussion we closely follow [50]. Let us introduce a set of subspaces of C n given 
by Vi = {ei,...,ej} C C n , where ej is the unit vector in the i th direction. A basis for 
the homology of the Grassmannian is determined by looking at planes C G G(k,n) whose 
intersection with our "reference" subspaces Vi is of a given specified dimension. The set of 
subspaces Vi is called a "flag" V = (Vi C V% C • • • C V n -i C V n = C n ). More precisely, we 
are interested in cycles of the form 

o- ai ,...,a k = {CE G(k, n) : dim(C n K-fc+i-J > i] (142) 

where {(ai, . . . , a&)} ranges over all nonincreasing sequences of integers between and n — k. 
These are called Schubert cycles. 

The intersection of Schubert cycles is well understood and has a fascinating connection 
to Young tableaux as one could have imagined from the restriction on the set of a/s. From 
the restriction it is clear that Schubert cycles are in one-to-one correspondence with Young 
tableaux. The whole intersection theory is encoded in the tensor product formula of irre- 
ducible representations of unitary groups. In other words, consider cycles and a\ where 
fj, and A are two Young tableaux. We know that 

H ® A = © p cJ iA p (143) 
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where the sum is over all Young tableaux and A are the Littlewood- Richardson coefficients. 
Back to Schubert calculus, the intersection of <r M and <j\ is given by the formula 

p 

Now we are ready to make the connection with the computation of residues we have to 
perform to get physical information. Note that when we set a given minor to zero we get a 
hypersurface in G(k,n). This hypersurface turns out to coincide with <j\ for some choice of 
flag^l The cohomology ring of the Grassmannian does not depend on the choice of the flag. 
Therefore when we talk about a given cycle we are implicitly taking about the corresponding 
class. 

The dimension of <j\ is one less than the dimension of the Grassmannian, i.e, of codimen- 
sion one. If we set another determinant to zero we get another representative of the class 
of <j\. Imposing both equations at the same time gives a cycles of codimension 2 which is 
homologous to the self-intersection of a\, i.e, o\. At this point we have to remember that 
we are not working directly on the full G(k,n). The constraints coming from momentum 
conservation restricts our computation to a smaller Grassmannian, i.e., G(k — 2, n — 4). This 
reduction was discussed in detail in section 12.21 

In order to compute residues we need to take enough self intersections of o\ so that we find 
cycles of dimension zero. Starting with o\ of dimension D — l where D — {k — 2) (n — k — 2) 
is the dimension of G(k — 2, n — 4), one has to compute erf. Recall that <j\ is associated 
with the tableaux consisting of a single box, i.e., the fundamental representation. Now we 
are interested in how many times the tableaux with k — 2 rows and n — k — 2 columns 
appearsj. In the case at hand the tableaux of interest corresponds to a cycle of codimension 
(k — 2) (n — k — 2), i.e, a point! The computation of the corresponding Littlewood- Richardson 
coefficient is a classic computation in enumerative geometry and gives 

# = U21...(fc-3)U* (145) 

W (n-Jfe-2)!(n-Jfe-l)!...(n-5)f 1 ' 

Let us test the formula against the cases we know. For k = 3 and all n we find # = 1. 
This is indeed the case as any choice of n — 5 determinants set to zero has only 1 solution 
since all equations are linear. For k = 4 and n = 8 we find # — 2. This is precisely the 
curious fact we noticed in previous sections. Also for k = 4 we have checked that when n = 9 
we get # = 5 solutions. The most impressive agreement is for the number of solutions in 
the N 3 MHV, i.e., k = 5, n = 10 amplitude. In this case, we have to find all solutions of a 
system of 9 cubic polynomial equations in 9 variables. With the help of Mathematica we 
have checked that there is precisely # = 42 solutions! 



5 We use the usual mathematical notation where cr a .o,...,o is denoted by a a . 

6 A simple way of computing the codimension of a cycle given its tableaux is to count the number of boxes. 
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9 Discussion 



We have identified an object that autonomously computes the leading singularities of quan- 
tum scattering amplitudes in an entirely different way than the usual local formulation of 
quantum field theory. There is a great deal left to be understood; we end with a discussion 
of open problems, speculations and directions for future work. 

9.1 Identifying Contours 

The most immediate open question we confront is how to identify contours associated with 
leading singularities for completely general k and n. As we have seen, the amplitudes are 
associated with equivalence classes of contours which are related by sums of residues that van- 
ish due to the global residue theorem. We have seen a beautiful pattern to these equivalence 
classes for NMHV amplitudes. Beginning with the 8 pt N 2 MHV amplitude the story gets 
even more interesting, involving "composite" residues in addition to the usual ones; in fact as 
we pointed out for large n = 2k, all the residues must be composite, so understanding their 
physical significance will surely be important. But obviously, apart from simply identifying 
the contours associated with known leading singularities, we need to find an independent 
principle which makes the correct associations. 

At 1-loop, all that is needed to completely determine the amplitude is a prescription for 
associating a contour T with every scalar box; our conjecture then says that the full amplitude 
is 



At higher loops, we need to learn both how to systematically associate residues with leading 
singularities as well as how to use the leading singularities to re-assemble the full amplitude. 
It is quite likely that these two steps are closely related to each other. A possible clue is that 
the procedure for identifying residues is highly reminiscent of the way in which the leading 
singularities themselves are defined. At 1-loop one begins with a loop integral, thought of 
as a contour integral over the whole real axis for the 4 loop momenta, and then defines the 
leading singularity as being associated with exactly the same integral but evaluated on the 
two T 4 's associated with the two solutions to the quadruple cut equations. At 2 loops, one 
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encounters "composite" leading singularities entirely analogous to our "composite" residues. 
This perhaps suggests that our expression for the leading singularities is in fact the reduction 
of some even larger integral representation of the full loop-level amplitude, and that the 
restriction to T u contours which picks out leading singularities also specifies the correct 
contours in our conjecture. There are many other indications of a beautiful structure in 
the full loop amplitudes. For instance at 1-loop, it is only the full amplitude (and not 
each individual scalar box) that exhibits the remarkable duality between Wilson loops and 
scattering amplitudes [Ml ESI [561 EH H31 EH EHl [601 [61] . 

9.2 Dual Superconformal Invar iance and Momentum Twistors 

As we have stressed our basic expression for the amplitude makes the action of the cyclic 
and parity symmetries manifest. As we have learned over the past few years, SYM am- 
plitudes also have a remarkable dual superconformal invariance, which together with the 
usual superconformal symmetry close into a Yangian algebra [62], demonstrating the same 
integrable structure in scattering amplitudes as has already been observed in the spectrum 
of anomalous dimensions [63] . While this symmetry has been understood as a consequence 
of a fermionic T-duality [59], it would be nice to understand it even more directly. Given 
that our proposal appears to unify all functions of the kinematical invariants that appear in 
scattering amplitudes in one object, it is very natural to explore whether it might make dual 
superconformal invariance "obvious". Indeed, perhaps the real interpretation of our object 
is as a generating function for all possible Yangian invariants which also knows about all the 
remarkable relations among them. 

Along these lines, it is perhaps worth pointing out some completely obvious re-writings 
of our initial twistor-space formula. We begin by writing the <5 4 ' 4, s as an integral over k 
auxiliary twistor variables Z a , and also writing each of the inverse determinant factors as a 
bosonic gaussian integral. This gives us a bosonic integral with a cubic action of the form 

C n , k (W a ) = J dX^dY^dZ a dC aa e* s , S = X^Y^C P:j+a ^ + W a C aa Z a (146) 

This has the form of an action for a spin chain, with interactions involving sets of k nearest- 
neighbors. 

Properly interpreting dual conformal transformations in our formulation should also help 
us understand how our picture is related to the recent work of Hodges [39]. His picture 
for NMHV amplitudes is of a volume integral in the dual twistor space, keeping the cyclic 
symmetries manifest, with the BCFW expression expressing the decomposition of this volume 
into tetrahedra, and unphysical poles associated with the internal faces of the polytope, 
which cancel in the sum representing the whole volume. This seems very closely related to 
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our picture in ordinary momentum or twistor space, where the amplitude is associated with 
a cyclically invariant homology class, even though any given contour integral representation 
does not have the symmetries manifest. On the other hand there are clear differences-our 
proposal also keeps parity manifest, and even for NMHV amplitudes, combines tree and loop 
information in the same object. Perhaps Hodges picture is more naturally associated with a 
"T-dual" of the twistor string theory. 

9.3 Pure Yang-Mills? 

One of the striking features of maximally supersymmetric theories, even at tree-level, is that 
all the different helicities are unified in a single supermultiplet and the amplitudes can be 
treated as a single object with external states depending on Grassmann parameters. In 
the usual formulation of scattering amplitudes, the relations amongst the different helicity 
amplitudes that follow from SUSY seem utterly mysterious without knowledge of their super- 
symmetric origin. However our conjecture, written in its non-supersymmetric form, retains 
the same essential feature provided by SUSY: the integrand is a universal object, and the 
different helicity configurations arise from performing the integration over different charts on 
the Grassmannian. This perhaps suggests that SUSY is not playing a particularly crucial 
role in the story, and that there may be an extension of our formula for computing the sub- 
leading singularities as well, corresponding at one-loop to the triangle, bubble and rational 
pieces. Finding a generating function for all rational terms in scattering amplitudes would be 
particularly exciting both from a purely theoretical point of view, since these objects clearly 
have a fascinating and deep structure, as well as for practical purposes, since computing the 
rational parts of 1-loop amplitudes is crucial for state-of-the-art NLO calculations in QCD 

9.4 Emergent Space-Time and Unitarity 

We are still missing a real understanding of the physics behind our conjecture. A clue is 
perhaps provided by the nature of the space in which it is formulated. Instead of space- 
time, for n-particle scattering the dual is naturally formulated in an n-dimensional space. 
Thinking along "holographic" lines, we may have expected a theory living on the boundary 
of spacetime or a close cousin like twistor space. However we are finding that our dual 
picture isn't associated with a space in which particles live at all; instead most of the action 
takes place in the Grassmannians G(k,n). All the kinematical information associated with 
spacetime goes into specifying special directions in this n-particle space, which pick out 
out a smaller G(k — 2,n — 4) Grassmannian naturally embedded in the larger one. Thus, 
while our duality is trivially "holographic" given that we are working directly with on-shell 



60 



variables, space-time emerges in a less directly "holographic" way than we are accustomed 
to in AdS/CFT, beginning instead from a picture where there is no space of any sort and 
the arena in which the dynamics takes place is determined by the number of particles. 

One might be tempted to loosely think of this n-particle space as analogous to a classical 
phase space, but the similarity seems superficial. Dynamics in classical phase space is local. 
By contrast, leading singularities are contour integrals in G(k — 2,n — 4), and they are 
associated with topological information. Indeed we can focus on the G(k — 2,n — 4) space, 
"punctured" by removing the zero locus defined by the vanishing of the minor factors; the 
contours defining the amplitudes are then associated with non-trivial homology classes in 
this space. This gives several equivalent representations to a given amplitude, which helps us 
see the emergence of space-time by enforcing the cancelation of unphysical poles and infrared 
equations at loop level. 

We emphasize the emergence of space-time since this is something one would expect 
to see in a dual theory of the S Matrix, but in fact it is possible that unitarity should 
also be thought of as emergent in this picture. After all, every discussion of scattering 
amplitudes begins by imagining all momenta as incoming, which immediately removes the 
unitary "matrix" interpretation of scattering! An even stronger indication that unitarity is 
not a fundamental ingredient in the dual we have been studying is that the object which is 
being directly computed is the leading singularity and not the unitarity cut. Consider the 
1-loop amplitude; given the leading singularities we can construct an object that has the 
correct quadruple cut, which (with maximal supersymmetry) also correctly reproduces all 
the subleading cuts. The first subleading cut is the triple cut, and the "double-cut" (which 
is the cut associated with unitarity) does not seem particularly privileged. The "double-cut" 
only acquires special significance working in (3, 1) signature with real momenta, since both 
the quadruple cut and triple cuts vanish in this case. 

A perhaps too-radical thought inspired by these observations is that the dual theory cares 
neither for space-time nor unitarity and computes "amplitudes" as determined by the leading 
singularities. The very special analytic structure of the "amplitudes" allows a local space- 
time description; restricting to spinor helicities satisfying a reality property A = A* gives 
(3, 1) signature and the possibility of a causal interpretation. The unitarity cut equations 
are satisfied as a rather incidental fact, and therefore the "amplitudes" have a probabilistic 
interpretation. In this way both space-time and the usual quantum mechanical rules may 
arise from a more primitive starting point. 
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9.5 M = 8 SUGRA 



Finally, we make some obvious comments about the extension of our work to M = 8 SUGRA 
[651 ESI EZ] , which is bound to be much more interesting [SHlEHlITOlinilTSlITallHIIHlEgEZl 
[781 [791 [801 EH [821 ESI IHH ESI EHl [87]. As stressed in [9], gravity amplitudes have a far richer 
structure than in Yang- Mills theory. They are governed by much larger "obvious" symmetries: 
instead of cyclic invariance, supergravity amplitudes have permutation symmetry, and instead 
of having the massless scattering amplitudes only defined at the origin of moduli space as 
for M = 4 SYM, the M = 8 SUGRA S Matrix is defined everywhere on moduli space and 
is non-trivially acted on by the £7(7) symmetry [91 [881 EH]. There are also indications of 
many more non-trivial relations between gravity amplitudes than exist for Yang-Mills. To 
begin with, gravity amplitudes are softer at complex infinity than gauge amplitudes, and this 
gives rise to extra relations for tree gravity amplitudes beyond those given by the BCFW 
recursions relations. A large number of such relations are needed in order to guarantee a 
remarkable property of loop-level gravity amplitudes: the soft limits are not renormalized 
beyond tree level. This is obvious from the usual local description but somewhat miraculous 
using on-shell methods, and its validity requires yet further non-trivial relations between tree 
gravity amplitudes. 

All of these observations suggest a far richer structure controlling gravity amplitudes 
already at tree level, and yet there has always been a note of discouragement in pursuing this 
thought: even the tree MHV gravity amplitudes look more complicated than the beautifully 
simple Parke- Taylor amplitude. We believe this is because we have yet to see supergravity 
amplitudes expressed in a form that makes the full permutation symmetry manifest, and that 
the analog of our M = 4 conjecture for M = 8 SUGRA will make it clear that supergravity is 
indeed the "simplest" theory in this dual formulation [9]. There are already very encouraging 
indications that the link representation for gravity gives rise to a remarkable new form for 
tree MHV gravity amplitudes [90J. 

There is an opportunity to study questions related to Af = 8 SUGRA by looking at 
subleading trace structures in J\f = 4 SYM. For instance, the soft X behavior of gravity 
amplitudes under BCFW deformation is also present for Af = 4 SYM amplitudes, when 
non-adjacent legs are deformed. All of the structure we have seen in this paper, including 
the remarkable identities between different BCFW representations of tree amplitudes, have 
instead held only for planar BCFW diagrams associated with deforming adjacent legs. This 
suggests that even Yang-Mills theory should exhibit a rich structure in its non-planar sector, 
and understanding this will be relevant to Af = 8 SUGRA as well. 
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A Relation to the Link Representation 



In this appendix we show how that our conjecture is related to what we called the "link 
representation" of scattering amplitudes in [18]. Amongst other things, this will allow us 
to interpret (and finally be rid of) one of the annoying features of representing scattering 
amplitudes in twistor space-the presence of "infinity twistor signs" [HI [191 EI] that obscure 
manifest conformal invar iance. It has always seemed vaguely like these signs are only present 
to ensure correct contour orientations, and we will see how this arises very precisely. 

Let us illustrate how this works for the simplest non-trivial case of the 6pt NMHV am- 
plitude. As explained in detail in [H], assigning i — 1, 3, 5 to W's and / = 2, 4, 6 to Z's, we 
can write the full superamplitude as M 6 = (1 + g 2 + g 4 )Us where the twistor diagram for Uq 
is 

H 3 




5 Gil (147) 

This also allows us to write an expression for the superamplitude in what we called "the link 
representation" in twistor space: 

U 6 = sgn((46)[13]) J d 9 c n U 6 ( Cli )e^ w ^ (148) 

with 



U(cn) = <5(c 25 ) 



1 



1 



1 



C45 c 65 C61 c 43( c 61 c 43 — C41CQ3) C21C23 



(149) 
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Since all the /tj,/// dependence of amplitudes in this representation is in the e lCuWxZl factor, 
it is trivial to go back to momentum space. Using the (say) ij representation for all particles 
we have 

U 6 = sgn((46)[13]) J d 9 c R U 6 (c n )5 2 (X t - C/,A/)5 2 (A/ + cnX^im + c u fu) (150) 

As explained in the text, the bosonic delta functions have linearized momentum conservation, 
since 

K - c R \i = 0, Aj + c u \i = XiXi + A/A/ = (151) 

Thus the 12 bosonic delta functions turn into the 4 delta functions for momentum conserva- 
tion plus 8 more delta functions restricting the c/«. Note that U% also contains a <5c 2 5 factor, 
so the remaining 8 en are completely determined by the 8 delta functions. For instance we 
can use the A 5 and A 2 equations to solve for c 65 , c 45 , c 2 i, c 23 , 

(54) _ (56) [23] _ [21] 

c 65 - py, c 45 - py, c 21 - c 23 - ^ (152) 

and then use, say, the A 4 , A 6 equations to solve for the rest of the c's, 

r (6|(p 5 +p 4 )|3] n (6|(p 5 +p 4 )|l] „ (4|(p 5 + P6)|3] „ (4|(p 5 +P6)|l] 
41 (46) [13] ' 43 (46) [31] ' 61 (64) [31] ' 63 (64) [13] 1 J 

From here it is straightforward to recover momentum space amplitudes. To obtain the 

(H 1 1 — ) amplitude, we integrate f d 4 fji, which simply gives 1, and set the remaining 

fji — > 0. Taking account the Jacobian in pulling out the overall momentum conserving 

5 function and inserting the solutions for the en into U 6 gives us the momentum space 

amplitude. We can also easily obtain the ( — I 1 h) amplitude, by integrating f d 4 fji, 

and setting fjj — > 0. The integrals yield a factor of det 4 (c/j), so we find 

tr+-+-+ - det 4 ( C /.)[/ + -+-+- - ( m+m y u +-+-+- 

m + 3)15)4 x „ ' n (154) 



[12] [23] (45) (56) (p 4 + p 5 + Pe) 2 (6|5 + 4|3] (4|5 + 6| 1] 
and of course 

M 6 -+-+-+ = (1 + g 2 + g')U % + - + - + - (155) 

As an aside we can quickly derive the identities that follow from the cyclic symmetries of the 
superamplitude. For instance cyclicity implies that M 6 H_+_+_ (Aj, \) = M & ~ + ~ + ~ + (X i+ i, X i+ i). 
Given the expressions we have just given for these two amplitudes we find a 6 term identity 

(46) 4 [13] 4 1 



[i + g' + g*) 



[12] [23] (45) (56) (p 4 + p 5 +p 6 ) 2 (6|5 + 4|3](4|5 + 6|1] 



[3|(2 + 4)|6) 4 1 



:i + 9 +9) [23] [34] (56) (61) (p 5 + p 6 + Pl f (1|6 + 5|4](5|6 + 1|2] (156) 
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Of course this is just one component of an identity as a function of the Grassmann parameters 
fji. This is precisely the same as the identity Mbcfw = -^p(bcfw) we mentioned in the 
introduction. 

At any rate, the whole tree amplitude is the sum of three terms; we would like to un- 
derstand them as all arising from a single object. The 6(025) is a clue; perhaps we could 
somehow replace this with a factor 1/(025) and think of U% as computing a residue on the 
pole where c 25 = 0. The sgn(46)[13] factor appears to be an annoying obstruction to this 
identification, but in fact as we will see it has an important role to play to ensure everything 
works out perfectly. 

We are actually very close to the desired result, but to see this, we need a small piece of 
inspiration: we have to recognize that the factor c 45 c 2 i is really the same as C45C21— C25C41 = c 63 
when evaluated on the support of 6(025), an d similarly for c 65 c 2 3 = c 65 c 2 3 — c 2 5C 63 = c 4 i, these 
are of the same form as the already appearing factor cuc^q — C16C34 = £25. In fact we were 
directly motivated to make this identification by studying the localization properties of this 
amplitude in twistor space, but it would take us too far afield to explain that here. For now 
we simply note that with the identification we can write Uq in a more symmetrical looking 
form as 

U 6 = [ d 9 c R e^ w ^ sgn((46)[13])^(c 25 ) 1 . j . (157) 

J C 4 iC 63 C2 5 C4iC 63 



where 



Cli — tijktlJKCjjCkK (158) 



This form already motivates the loose association 

1 



Uq ~ Res 



(159) 

J C 25 =0 



_C25C6lC4iC25C63C4i_ 

which we will make completely precise in a moment. Note however before proceeding that 
the factors in the denominator are in perfect agreement with the consecutive minors of of the 
matrix 

(c 2 i 1 c 23 c 25 \ 
C41 c 43 1 C45 . (160) 
cei c 63 c 65 1 / 

Namely 

C41 = (612), c 63 = (234), c 25 = (456), c 25 = (123), c 41 = (345), c 63 = (561) (161) 

In order to do complete the transformation of our integral involving 6(025) into the calculation 
of a residue, let us go back to momentum space, and repeat the steps we described in section 
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2 in the analysis of our formula. Looking at the space of link variables on the support of the 
bosonic delta functions, 

Ai - cnh = 0, \i + cnXi = (162) 

we have 12 equations, but 4 are redundant given 4-momentum conservation, so there are 
8 equations for 9 unknowns. The set of en solving these equations lie on a line in the 9 
dimensional cu space. We can parametrize these solutions e.g. by solving for all the c's as a 
function of one of them, or we can parametrize them more elegantly as 

ch(t) = c*j + €i jk e IJK \jk]{JK)T (163) 

where c*j is any specific solution. Here we use the fact that for three two-dimensional vectors 
e.g. Xi,j,k, 

€ U kXi{JK) = (164) 

We can then write the product of the S 2 functions as the momentum conserving delta function 
multipled by an integral that localizes the cu to this line: 

\\?>\\-cn\i) Y[6 2 (\i + cn\i) ^ 4 £U) f dr5 9 (c H - c«(r)) (165) 

i I a J 

(it is easy to check that the Jacobian is 1). 

Using this we can re- write Ug~ + ~ + ~ , explicitly factoring out the momentum conserving 
delta function, as 

Ut ,- + - = f dT fW^'l ( 166) 

Now, something very nice happens, which allows us to finally interpret-and be rid of!-the 
pesky sign factors! We can always shift the origin of r so that 

c 25 (r) = (46)[13]r (167) 

Then doing the r integral simply gives us 

U?- + - + - = 5\Y Pa ) - \ 1 x — (168) 

6 ^^(46) [13] (c 4 ic 6 3C 2 5C4ic 63 )(r = 0) V ; 

But we can immediately recognize 



1 



(46) [13] 4 (rH c 25 (r) 



(169) 
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and therefore we have finally identified Uq as a residue, which 

Ut + - + - = ^ 4 (E^) / 7 ~ ~ - v S ( 17 °) 

^ Jc 25 (t)=0 (C25C63C4lC2 5 C 6 3C4l)(r) 

which we can re-write in precisely the form of our conjecture using equation (I16ip as 

U " ++ = 64 % Pa) £ {T)=0 [(123)(234)(345)(456)(561)(612)](r) 

Note that in the end the role of the annoying sgn factors was just to tell us to consistently 
define all residues with the same orientation! 

As discussed in the text, the BCFW form of the full tree amplitude is given by a contour 
integral enclosing the poles where (612), (234), (456) vanish, while the P(BCFW) form of 
the amplitude used the other three poles, and the remarkable 6 term identity equating these 
two forms just follows from Cauchy's theorem. 

B An Explicit 7 particle Computation 

In this appendix we give an explicit derivation of a class residues for the 7 particle NMHV am- 
plitude, associated with the three-mass boxes. This only involves solving linear equations and 
is completely straightforward. We include it here so the reader can transparently see what is 
involved; as mentioned in the text for the cases with higher k and n we leave the computation 
of residues to Mathematica. Up to cyclic permutations, there is a unique three-mass box for 
the 7 particle amplitude. We will compute the box with particles (2), (3,4), (5,6), (7, 1) at 
each of the four corners. We find it convenient to use the gauge fixing 

/ 1 c u ci 5 cie ci 7 \ 
C = 1 c 24 c 25 c 26 c 27 (172) 
\ 1 c 34 c 35 c 36 c 37 / 

We now have to examine the linear equations Aj — cuXj = 0, Xj + CjjAj = 0, as well as 
the linear equations involved with setting minors to zero. Here I = 1, 2, 3 and i = 4,5, 6, 7. 
We can solve these equations in any order we wish, and for convenience we will first impose 
the A equations, then the equation setting the minors to zero, and finally the A equation. 

Let us begin with the Aj equations: 



Aj — CjiAi + q 2 A 2 + q 3 A 3 (173) 

These are two equations for three unknowns. We parameterize the solution as 

(i 3) - (2 3)r; , . . . (i 1) - (2 l)n , . 

Ql(Ti) (13) ' °afa)=n> <tafo)= \ L-i- 1 \ (174) 
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The three-mass box coefficient will correspond to the residue that occurs when the determi- 
nants (234) and (712) vanish. This allows us to solve for 

T4= (23)' T7= (21) (175) 

Next let's impose the \j + cu\ = equation. This constraint translates into two equa- 
tions which fix the remaining two r's and introduce a dependence on the anti-holomorphic 
brackets-one might have expected six equations, but four turn into momentum conservation 
on the support of the other two. For example we can choose to impose 

A 2 = r <**> ( 176 ) 

#{1,2,3} 

and since A are two-component objects these constitute two equations. The solutions are 

[76](71) [64](43) [62] 

5 " [65] (21) [65] (23) [65] 1 ' 

[75] (71) [54] (43) [52] 

6 ' [56] (21) [56] (23) rccl ' ' ' 



To obtain our final expression for the box coefficient, we need to evaluate the determinants 
(2,3,4), (3,4,5), (4,5,6), and (5,6,7) on the pole, and compute the associated Jacobians 
from the delta function constraints, the poles, and the transition to the standard momentum 
conservation delta function. Plugging in for Tj in the determinants, we obtain 

(43)[6|1 + 7|2) 

(3 ' 4 ' 5) " [65](32)(21) (179) 

f4 5 6 x (21)([2|5 + 6|4)(32) + tjf'(43)) + [7|5 + 6|4)(32)(71) 

1 ' ' ' [65](21)(31)(32) 1 1 

( ,, 7] (21)(-[2|5 + 6|7)(32) - [4|5 + 6|7)(43)) - tj, 3l (32)(71) 

(5 ' 6 ' 7) " [65] (21) (31) (32) (181) 

(6,7,1) = ™ + 4 ' 2 > m 
V ' ' ' [65] (32) (21) v ; 

while (1, 2, 3) = 1 and (2, 3, 4) = (7, 1, 2) = 0. A short computation shows that the Jacobian 
takes the very simple form 

J = [56] (12) (23) (183) 
Thus, our final expression for the three-mass box coefficient at n = 7 is 

L 3J - 5 I ^ * ) (3,4,5)(4,5,6)(5,6,7)(6,7,1) 



68 



In order to extract the box coefficient corresponding to various helicity assignments for the 
external particles, we simply perform the appropriate grassman integrations. For instance, if 
want the three-mass box coefficient for the split helicity configuration M(l~2~3~4 + 5 + 6 + 7 + ), 
then because of our particular gauge fixing choice, we can simply set all the 77's to zero, giving 

LW} " = 64 ^ ) (3,4,5)(4,5,6K5,6,7)(6,7,1) (185) 



i=l 



If instead we want the three-mass box coefficient for the mostly alternating configuration 
M(l~2 + 3 _ 4 + 5~6 + 7 + ), then we need to flip the helicities of particles 2 and 5 relative to the 
split helicity case. For this reason we grassman integrate L37 over J d 4l fj2d 4l fj 5 e' ll2fl2+V5fl5 , and 
then set all the 77's and rfs to zero. Because of the simple dependence on the grassman 
variables, this integral simply brings down a factor of c 25 , so 



1 



L{im,7 - 5 4 (E A ^) (3j 4j 5) (4,5,6X5,6,7) (6, 7,1) (186) 



,i=i 



[6|7 + 2|1)(23) + [6|4|3)(21) 
25 " [65] (12) (23) 1 ' 

where we have solved for C25 in terms of the Tj. Box coefficients for all other helicity combi- 
nations of the NMHV seven-particle can be computed in this fashion. 



C IR equations at 1-loop 

Here we give a short description of IR equations at 1-loop, for more details see [HI [53] . 
One-loop partial amplitudes in Yang Mills are IR divergent. Due to color ordering, the 
divergences can only depend on the momentum of two consecutive external particles and 
it must be proportional to the tree-level amplitude. More explicitly and in dimensional 
regularization, 

1 n 

<- loop | IR = -j % E(- s vm) e Mr e . (188) 
i=i 

As discussed in previous sections, one-loop amplitudes can be written as linear combinations 
of scalar box integrals. Each integral can have IR divergences that depend on the kinematical 
invariants of the particular box under consideration. Only the four-mass box integral is 
completely finite. 

Let us denote each box integral by I(Kx, K 2 , K 3 , K 4 ) as in section 16.11 and the main 
kinematical invariants by s = {Ki + K 2 ) 2 and t = (K 2 + K 3 ) 2 . Each Ki is the sum of 
consecutive external momenta. If Ki is equal to the momentum of a single particle then it is 
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null, if it is the sum of two or more then we say that it is massive. We denote null momenta 
by lower case letters (p,q,r) while massive momenta by capital letters (P,Q,R). Using 
our notation one-mass, two-mass-easy, two-mass-hard and three-mass boxes are respectively 
given by: I(p, q, r, P), I(p, P, q, Q), I(p, q, P, Q), and I(p, P, R, Q). 

We now list the IR divergent structure of each of the four classes of IR singular boxes: 

I(p,q,r,P) = -I((- s )- + (-t)--(-P 2 )-), 
I(p,P,q,Q) = -^((- S )- e + (-t)- e -(-P 2 )- e -(-g 2 )- e ), 
I(p,q,P,Q) = -^Q(-3)- e + H)- e -i(-P 2 )- £ -i(-g 2 )- e 

I( P ,P,R,Q) = -I Q(- a )- + -i(-^)--i(-g»)-^. (189) 

We would like to propose a distinction between two classes of IR equations. There is 
the class that relates the one-loop coefficients to the tree-level amplitude and the class that 
relates one-loop coefficients purely among themselves. In the following we will argue that 
this separation is meaningful. Moreover, we will find that the latter class is a subset of more 
purely one-loop identities which turn our to be derivable from a generalization of the residue 
theorem. 

The first class is obtained by collecting all terms in the expansion in scalar box integrals 
that are proportional to a given — \{— Si t i+i) e divergence. The sum of the corresponding 
coefficients weighted by the numerical factors in (11891) must then be equal to M* ree . The 
second class is obtained by collecting all terms in the expansion in scalar box integrals that 
are proportional to — \(— (pi+pi + \ + . . .+p i+m ) 2 ) e with m > 1. The sum of the corresponding 
coefficients weighted by the numerical factors in H189I) must now vanish. 

The first class is used to obtain the |[BCFW + P(BCFW)] form of recursion relations for 
tree amplitudes. An example of the second class which we use in the text is provided by the 
seven particle amplitude. Consider for instance the equation coming from £123, 

= -P(l)( 2 )(3)(4567) + 5(7) (1) (23) (456) + #(3) (4) (567) (12) ~ |-B(4)(5)(67)(123) ~ | #(6) (7) (123) (45) 
~-B(4)(56)(7)(123) + #(1) (23) (4) (567) + #(3) (456) (7) (12) + § #(3) (45) (67) (12) + §-B(l)(23)(45)(67) • 

As we discuss in the text, this highly non-trivial identity is exactly equivalent to a residue 
theorem. 
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